AD-A 


David  Taylor  Research  Center 


Bethesda,  MD  20084-5000 


DTRC-90/009  April  1990 

Propulsion  and  Auxiliary  Systems  Department 
Research  and  Development  Report 


one  FILE  COPY 


Propulsive  Efficiencies  of  Magnetohydrodynamic 
Submerged  Vehicular  Propulsors 

by 

S.  H.  Brown,  J.S.  Walker,  N.A.  Sondergaard, 

P.J.  Reilly,  and  D.E.  Bagley 


DTIC 

o^YoTsn 


Approved  tor  public  release;  distribution  is  unlimited. 


CODE  Oil  DIRECTOR  OF  TECHNOLOGY,  PLANS  AND  ASSESSMENT 
12  SHIP  SYSTEMS  INTEGRATION  DEPARTMENT 

14  SHIP  ELECTROMAGNETIC  SIGNATURES  DEPARTMENT 

15  SHIP  HYDROMECHANICS  DEPARTMENT 

16  AVIATION  DEPARTMENT 

17  SHIP  STRUCTURES  AND  PROTECTION  DEPARTMENT 

18  COMPUTATION,  MATHEMATICS  &  LOGISTICS  DEPARTMENT 

19  SHIP  ACOUSTICS  DEPARTMENT 

27  PROPULSION  AND  AUXILIARY  SYSTEMS  DEPARTMENT 

28  SHIP  MATERIALS  ENGINEERING  DEPARTMENT 


DTRC  ISSUES  THREE  TYPES  OF  REPORTS: 

1 .  DTRC  reports,  a  formal  series,  contain  information  of  permanent  technical  value. 
They  carry  a  consecutive  numerical  identification  regardless  of  their  classification  or  the 
originating  department. 

2.  Departmental  reports,  a  semiformal  series,  contain  information  of  a  preliminary, 
temporary,  or  proprietary  nature  or  of  limited  interest  or  significance.  They  carry  a 
departmental  alphanumerical  identification. 

3.  Technical  memoranda,  an  informal  series,  contain  technical  documentation  of 
limited  use  and  interest.  They  are  primarily  working  papers  intended  for  internal  use.  They 
carry  an  identifying  number  which  indicates  their  type  and  the  numerical  code  of  the 
originating  department.  Any  distribution  outside  DTRC  must  be  approved  by  the  head  of 
the  originating  department  on  a  case -by -case  basis. 


NOW -DTNSROC  5602/51  (Rev  2-88) 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 


1*.  REPORT  SECURITY  CLASSIFICATION 

Unclassified 


lb.  RESTRICTIVE  MARKINGS 


24.  SECURITY  CLASSIFICATION  AUTHORITY 


3.  DISTRIBUTION/AVAILABILITY  OF  REPORT 


2b.  DECLASSIFICATION/DOWNGRADING  SCHEDULE 


Approved  for  public  release;  distribution  is  unlimited. 


4.  PERFORMING  ORGANIZATION  REPORT  HUMBERTS) 

DTRC-90/009 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6*.  NAME  OF  PERFORMING  ORGANIZATION 

David  Taylor  Research  Center 


6b.  OFFICE  SYMBOL 
(If  applicabla) 

Code  271 


7*.  name  of  monitoring  organization 


6c.  ADDRESS  ( City,  Sam,  and  ZIP  Coda) 


Bethesda,  Maryland  20G84-5000 


7b.  ADDRESS  (City,  Saw,  and  ZIP  Coda) 


64.  NAME  OF  FUNDING/SPONSORING 
ORGANIZATION 

(See  reverse  side) 


6b.  OFFICE  SYMBOL 
(If  applicabla) 

DTRC0113 


6.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


ADDRESS  (City.  Stata.  and  ZIP  Coda) 

DTRC 

Bethesda,  Maryland  20034-5000 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

ELEMENT  NO. 

NO. 

NO. 

ACCESSION  NO. 

62936N 

1-2712-131 

ZF664 12001 

500540 

11.  TITLE  {Induct*  Security  CiM$sJtlcaticn) 

Propulsive  Efficiencies  of  Magnetohdrodynamic  Submerged  Vehicular  Propulsors 


12.  PERSONAL  AUTHOFKS) 

Samuel  H.  Brown,  John  S.  Walker,  Neal  A.  Sondergaard,  Patrick  J.  Reilly,  and  David  E.  Bagley 


134.  TYPE  OF  REPORT 

Final 


13b.  TIME  COVERED 
FROM _ TO 


14.  DATE  OF  REPORT  ( YEAH  UONTH.  DAY) 

April  1990 


15.  PAGE  COUNT 

64 


i6.  supplementary  NOTATION  -plus  WOrk  was  a  cooperative  effort  with  the  Depannent  of  Maechanical  and  Industrial  Engineering 
at  the  University  of  Illinois,  Urbana  HI.  61801 _ _ 


17.  COSATl  CODES  | 

FIELD 

GROUP 

SUB-GROUP 

16  SUBJECT  TERMS  (Cordnua  on  ravaraa  *  naoaaaaiy  and  kfandty  by  block  nurnbaf) 


Magnetohydrodynamic  propulsion,  marine  propulsion,  seawater  pump, 
superconducting  magnets 


16.  ABSTRACT  (Condnua  on  ravaraa  d  naoaaaary  and  kiandty  by  btock  ntMitbat) 

4- 

Magnetohydrodynamic  (MHD)  ship  propulsion  is  the  process  of  propelling  a  vehicular  structure  by  a  seawater  elec¬ 
tromagnetic  pump.  This  propulsion  system  can  be  applied  to  a  surface  ship  or  a  submerged  vehicle;  however,  in  this 
work  only  submerged  vehicles  at  depths  where  wave  effects  can  be  neglected  were  considered.  Although  a  number  of 
different  arrangements  for  a  MHD  propulsion  system  are  possible,  the  general  characteristics  of  such  systems  are  most 
easily  determined  by  a  simple,  ideal  MHD  rectangular  duct  of  constant  cross-sectional  area.  A  mathematical  model  was 
developed  at  the  David  Taylor  Research  -Center  (DTR€ifrfor  calculating  the  propulsive  efficiencies  of  such  a  rectangular 
duct  propelling  a  submerged  vehicle.  Numerical  propulsive  efficiencies  are  presented  in  terms  of  many  different  param¬ 
eters.  Assumptions  were  generally  made  in  the  model  that  tend  to  maximize  the  propulsive  efficiency  of  the  'MHD  sys¬ 
tem.  Thus,  the  propulsive  efficiencies  calculated  from  the  model  overestimate  the  efficiencies  of  the) corresponding  real 
MHD  propulsion  system.  These  numerical  results  can  be  used  for  engineering  estimates  of  the  propulsive  efficiencies  of 
real  MHD  propulsion  systems. 

_ (Continued  on  reverse  side) 


20.  OISTFUBUnON/AVAILABIUTY  OF  ABSTRACT 

□  UNCLA8SIFIE0AJNLIMTTED  IS  SAME  AS  RPT  □  OTIC  USERS 


21 ,  ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 


224.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Samuel  H  Brown 


22b.  TELEPHONE  (fnduda  Araa  Coda) 

(301  >-267-2146 


22c.  OFFICE  SYMBOL 

Code  2712 


DO  FORM  1473,  JUN  86 


Previous  editions  are  obsolete. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


01 02-LF-0 14-6602 


UNCLASSIFIED 


Block  8a: 

DTRC  Independent  Research  Program,  Director  of  Naval  Research; 
Submarine  Block  Program,  Office  of  Naval  Technology 


Block  13  (Continued): 

A  comparison  of  the  ideal  model  MHD  rectangular  duct  propulsive  efficiencies  with  the  corresponding  effi¬ 
ciency  of  a  conventional  submerged  vehicle  propeller  (order  of  0.85)  generally  shows  the  MHD  efficiencies  to  be 
less  under  most  conditions.  The  cases  where  MHD  propulsive  efficiencies  are  in  ranges  that  approach  that  of  a  pro¬ 
peller  correspond  to  MHD  duct  volumes  greater  ihen  300  m3,  with  small  surface-area-to-volume  ratios,  or  the  mag¬ 
netic  inductions  are  in  the  order  of  20  tesla  or  more.  Magnetic  inductions  of  about  6  to  10  tesla  in  large  volumes  of 
space,  such  as  MHD  ducts,  are  probably  about  the  highest  field  strengths  possible  with  the  use  of  present-day  super¬ 
conducting  magnet  construction  technology;  to  produce  even  these  fields  will  require  significant  engineering  effort 
in  magnet  construction.  Therefore  an  important  area  of  concern  for  MHD  propulsion  is  the  construction  of  practical 
superconducting  magnets  and  support  structures  of  reasonable  mass. 

According  to  the  calculations  using  the  mathematical  model,  a  weakness  of  MHD  propulsion,  for  a  viable  sub¬ 
merged  vehicle  propulsion  system,  is  the  generally  low  propulsive  efficiencies  of  the  system  at  duct  magnetic  induc¬ 
tions  of  10  tesla  or  less,  when  compared  with  conventional  propellers.  Furthermore,  the  general  trend  is  for  propul¬ 
sive  efficiencies  to  decrease  as  vehicle  speed  increases.  Future  work  will  address  the  degrading  effects  on  the  rectan¬ 
gular  MHD  channel,  caused  by  flinging  magnetic  induction  and  current  distributions  at  the  duct  ends. 

A  very  important  and  potentially  promising  aspect  of  MHD  propulsion  is  the  possible  reduction  of  the  sub¬ 
merged  vehicle’s  acoustic  noise  signature.  Unlike  a  propeller,  the  MHD  duct  has  no  mechanical  moving  parts.  This 
may  make  the  MHD  propulsor  intrinsically  quieter  (acoustically)  than  a  rotating  propulsor.  Also,  the  replacement  of 
the  propeller  by  a  MHD  thruster  eliminates  the  requirement  for  a  high  pressure  rotating  shaft  seal.  Therefore,  in 
spite  of  the  model’s  lower  predicted  propulsive  efficiencies,  the  acoustic  aspects  of  MHD  propulsive  should  be  in¬ 
vestigated  in  detail,  and  the  possibility  of  designing  and  constructing  a  practical,  acoustically  quiet  submerged  ve¬ 
hicle,  even  at  the  possible  expense  of  propulsive  efficiency,  should  be  explored. 


Accession  For 

NTIS  GRA&I 

DTIC  TAB 

□ 

Unannounced 

□ 

Justlf lcatl  on — 

By - - 

- - 

Distribution/  _ 

Availability  Codes_ 
(Avail  and/or 
Special 


\ 

I 


CONTENTS 


Page 


Abstract .  1 

Administrative  Information  .  2 

Introduction  .  2 

History  .  2 

MHD  Propulsion  Concept  . 2 

MHD  Project  Goals  .  2 

Goals  of  this  Technical  Report .  3 

Traditional  Propeller  Theory  and  Approximate 

Extension  to  MHD  Ducts .  4 

Background .  4 

Approximate  Relationship  of  Propeller  to  MHD  Duct .  10 

Mathematical  Model  of  Magnetohydrodynamic  (MHD)  Ship 

Propulsion  with  Propulsive  Efficiency .  13 

Fundamental  Equations  of  MHD  Propulsion .  13 

Conservation  of  Linear  Momentum .  14 

Principle  of  Action  and  Reaction .  16 

Rate  of  Change  of  Kinetic  Energy  (Energy  Theorem) .  17 

Electromagnetic  Energy  Input  .  19 

Efficiencies  for  MHD  Propulsion .  20 

Numerical  Results .  22 

Discussion  and  Conclusions .  51 

Acknowledgments .  53 

Appendix  A.  Three-Dimensional  Effects  .  54 

References .  57 

FIGURES 

1 .  Stream  tube  for  propeller  system .  5 

2.  Control  volume  in  space .  6 

3.  Control  volume  for  propeller  mathematical  model .  7 

4.  Simple  MHD  rectangular  duct .  10 

5.  MHD  rectangular  duct .  13 

6.  Propulsive  magnetohydrodynamic  rectangular  duct  configuration .  23 

7.  Power  consumption  versus  submerged  vehicle  speed .  24 


iii 


DTRC-90/009 


8.  Propulsive  efficiency  of  square  MUD  duct  versus  vehicle  speed 

in  meters  per  second .  26 

9.  Thrust  required  to  propel  submerged  vehicle  at  constant  speed  V .  28 

10.  Propulsive  efficiency  of  square  MHD  duct  versus  magnetic 

induction  in  tesla .  29 

1 1 .  Propulsive  efficiency  versus  electrical  conductivity  of  seawater .  30 

12.  Propulsive  efficiency  of  square  MHD  duct  versus  duct  cross-sectional  area .  31 

13.  Propulsive  efficiency  of  square  duct  versus  duct  length .  32 

14.  Propulsive  efficiency  T|i  versus  duct  height  (duct  width  w  =  1  m) .  34 

15.  Propulsive  efficiency  T)i  versus  duct  height  (duct  width  w  =  2  m) .  35 

16.  Propulsive  efficiency  rj  j  versus  duct  height  (duct  width  w  -  4  m) .  36 

17.  Propulsive  efficiency  T|i  versus  duct  height  (duct  width  h>  =  8  m) .  37 

18.  Propulsive  efficiency  T|i  versus  duct  height  (duct  width  w  =  16  m) .  38 

19.  Streamlines  due  to  uniform  body  force  JB  distributed  over 

rectangular  volume  with  no  duct  walls  .  55 

TABLES 

1 .  Baseline  numerical  results  for  square  rectangular  duct .  25 

2.  Propulsive  efficiency  T^xlO2  versus  submerged  vehicle  speed  and  length 

of  square  duct  at  2  tesla .  39 

3.  Propulsive  efficiency  TijXlO2  versus  submerged  vehicle  speed  and  length 

of  square  duct  at  6  tesla .  40 

4.  Propulsive  efficiency  TjiXlO2  versus  submerged  vehicle  speed  and  length 

of  square  duct  at  10  tesla .  41 

5.  Propulsive  efficiency  TjiXlO2  versus  submerged  vehicle  speed  and  length 

of  square  duct  at  6  tesla  .  42 

6.  Propulsive  efficiency  T^xlO2  versus  submerged  vehicle  speed  and  length 

of  square  duct  at  10  tesla .  43 

7.  Propulsive  efficiency  TfiXlO2  versus  electrical  conductivity  and 

submerged  vehicle  speed .  44 

8.  Propulsive  efficiency  T^xlO2  versus  electrical  conductivity  and 

magnetic  induction  of  square  duct .  45 

9.  Propulsive  efficiency  TiiXlO2  versus  hydrodynamic  resistance  of 

submerged  vehicle  (hydrodynamic  resistance  R  =  1000  V2  newtons) .  46 

10.  Propulsive  efficiency  TjiXlO2  versus  hydrodynamic  resistance  of 

submerged  vehicle  (hydrodynamic  resistance  R  =  2000  V 2  newtons) .  47 


IV 


DTRC-90/009 


11.  Propulsive  efficiency  T|iXl02  versus  hydrodynamic  resistance 

of  submerged  vehicle  (hydrodynamic  resistance  R  =  3000  V 2  newtons) 

12.  Propulsive  efficiency  ffiXlO2  versus  hydrodynamic  resistance 

of  submerged  vehicle  (hydrodynamic  resistance  R  =  4000  V2  newtons) 

13.  Propulsive  efficiency  T)iXl02  versus  hydrodynamic  resistance 

of  submerged  vehicle  (hydrodynamic  resistance  R  =  5000  V 2  newtons) 


DTRC-90/009 


V 


ABSTRACT 


Magnetohydrodynamic  (MHD)  ship  propulsion  is  the  process  of  propelling  a  vehic¬ 
ular  structure  by  a  seawater  electromagnetic  pump.  This  propulsion  system  can  be 
applied  to  a  surface  ship  or  a  submerged  vehicle;  however,  in  this  work  only  submerged 
vehicles  at  depths  where  wave  effects  can  be  neglected  were  considered.  Although  a 
number  of  different  arrangements  for  a  MHD  propulsion  system  are  possible,  the  general 
characteristics  of  such  systems  are  most  easily  determined  by  a  simple,  ideal  MHD  rect¬ 
angular  duct  of  constant  cross-sectional  area.  A  mathematical  model  was  developed  at 
the  David  Taylor  Research  Center  (DTRC)  for  calculating  the  propulsive  efficiencies  of 
such  a  rectangular  duct  propelling  a  submerged  vehicle.  Numerical  propulsive  efficien¬ 
cies  are  presented  in  terms  of  many  different  parameters.  Assumptions  were  generally 
made  in  the  model  that  tend  to  maximize  the  propulsive  efficiency  of  the  MHD  sys¬ 
tem.  Thus,  the  propulsive  efficiencies  calculated  from  the  model  overestimate  the 
efficiencies  of  the  corresponding  real  MHD  propulsion  system.  These  numerical  results 
can  be  used  for  engineering  estimates  of  the  propulsive  efficiencies  of  real  MHD  propul¬ 
sion  systems. 

A  comparison  of  the  ideal  model  MHD  rectangular  duct  propulsive  efficiencies  with 
the  corresponding  efficiency  of  a  conventional  submerged  vehicle  propeller  ( order  of 
0.85)  generally  shows  the  MHD  efficiencies  to  be  less  under  most  conditions.  The  cases 
where  MHD  propulsive  efficiencies  are  in  ranges  that  approach  that  of  a  propeller  corre¬ 
spond  to  MHD  duct  volumes  greater  then  300  m3,  with  small  surface-area-to-volume 
ratios,  or  the  magnetic  inductions  are  in  the  order  of  20  tesla  or  more.  Magnetic  induc¬ 
tions  of  about  6  to  10  tesla  in  large  volumes  of  space,  such  as  MHD  ducts,  are  probably 
about  the  highest  field  strengths  possible  with  the  use  of  present-day  superconducting 
magnet  construction  technology;  to  produce  even  these  fields  will  require  significant  en¬ 
gineering  effort  in  magnet  construction.  Therefore  an  important  area  of  concern  for 
MHD  propulsion  is  the  construction  of  practical  superconducting  magnets  and  support 
structures  of  reasonable  mass. 

According  to  the  calculations  using  the  mathematical  model,  a  weakness  of  MHD 
propulsion,  for  a  viable  submerged  vehicle  propulsion  system,  is  the  generally  low  pro¬ 
pulsive  efficiencies  of  the  system  at  duct  magnetic  inductions  of  10  tesla  or  less,  when 
compared  with  conventional  propellers.  Furthermore,  the  general  trend  is  for  propulsive 
efficiencies  to  decrease  as  vehicle  speed  increases.  Future  work  will  address  the  degrad¬ 
ing  effects  on  the  rectangular  MHD  channel,  caused  by  fringing  magnetic  induction  and 
current  distributions  at  the  duct  ends. 

A  very  important  and  potentially  promising  aspect  of  MHD  propulsion  is  the  possi¬ 
ble  reduction  of  the  submerged  vehicle's  acoustic  noise  signature.  Unlike  a  propeller,  the 
MHD  duct  has  no  mechanical  moving  parts.  This  may  make  the  MHD  propulsor  intrinsi¬ 
cally  quieter  (acoustically)  than  a  rotating  propulsor.  Also,  the  replacement  of  the 
propeller  by  a  MHD  thruster  eliminates  the  requirement  for  a  high  pressure  rotating  shaft 
seal.  Therefore,  in  spite  of  the  model's  lower  predicted  propulsive  efficiencies,  the  acous¬ 
tic  aspects  of  MHD  propulsive  should  be  investigated  in  detail,  and  the  possibility  of 
designing  and  constructing  a  practical,  acoustically  quiet  submerged  vehicle,  even  at  the 
possible  expense  of  propulsive  efficiency,  should  be  explored. 
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INTRODUCTION 


HISTORY 

Because  nuclear  reactors  use  circulating  liquid  metals  as  a  heat  transfer  medium,  a 
substantial  effort  in  the  past  has  been  devoted  to  developing  liquid  metal  electromagnetic 
pumps.1-3  These  pumps,  besides  being  quite  reliable,  are  unusual  in  their  design,  since 
they  have  essentially  no  moving  parts  but  the  electrical  conducting  fluid  being 
pumped.  The  high  electrical  conductivity  of  the  liquid  metals  permits  pump  power  effi¬ 
ciencies  of  nearly  20%  under  certain  conditions. 

The  direct  current,  Faraday-type  pump  is  the  most  elementary  electromagnetic 
pump.  The  liquid  metal  is  pumped  through  a  duct  of  rectangular  cross  section.  Two  op¬ 
posite  walls  of  the  duct  serve  as  electrodes  which,  when  maintained  at  an  electric 
potential,  cause  a  high  current  to  flow  through  the  fluid.  A  strong  magnetic  induction  is 
passed  through  the  other  two  walls.  Thus,  the  current  direction  and  the  field  are  mutually 
perpendicular  to  each  other  and  perpendicular  to  the  fluid  flow  in  the  pump  caused  by  the 
Lorentz  body  force. 

MHD  PROPULSION  CONCEPT 

In  the  mid  1960’s  the  electromagnetic  pump  concept  was  converted  into  an  overall 
marine  propulsion  concept  by  Stewart  Way.4-7  The  basic  idea  was  to  pass  a  current 
through  seawater  in  the  presence  of  a  magnetic  induction,  in  order  to  accelerate  the  sea¬ 
water  and  to  provide  a  reaction  force  on  the  marine  vehicle  holding  the  magnet  and 
electrodes.  To  our  knowledge  this  magnetohydrodynamic  (MHD)  propulsion  concept  has 
not  changed  significantly  since  then.  The  concept  can  be  applied  to  surface  ships  and 
submerged  marine  vehicles.  A  number  of  references  are  available  in  the  literature  on 
MHD  propulsion. 8-17 

MHD  PROJECT  GOALS 

The  objectives  of  the  MHD  marine  propulsion  research  and  development  project 
initiated  at  the  David  Taylor  Research  Center  (DTRC)  were  to  examine  fundamental  prin¬ 
ciples  underlying  these  MHD  propulsion  systems  and  then  to  produce  critical  scientific 
and  engineering  data  and  conceptual  design  information.  This  program  is  being  coordi¬ 
nated  with  the  Defense  Advanced  Research  Projects  Agency  (DARPA)  MHD  Propulsion 
Program,  a  cooperative  effort  with  the  Naval  laboratories,  Argonne  National  Laboratory, 
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industry,  and  academia.  All  the  parameters  that  characterize  MHD  propulsion  systems 
are  to  be  extensively  studied.  This  approach  will  establish  a  baseline  group  of  exper¬ 
imental  data  and  theoretical  numerical  results  to  assess  the  feasibility  of  marine  MHD 
propulsion  for  Naval  surface  ships  or  submerged  vehicles. 

It  is  well  known  that,  in  general,  MHD  propulsion  systems  have  significantly  lower 
propulsive  efficiencies  than  conventional  propellers  under  basically  the  same  propulsion 
conditions.  However,  it  is  generally  assumed  that  these  systems  will  generate  much  less 
acoustic  noise  than  propellers,  thus  making  MHD  systems  more  difficult  to  detect  acous¬ 
tically  than  conventional  propulsion  systems.  Along  with  the  elimination  of  rotating 
seals,  gears  and  shafting,  this  noise  reduction  is  the  major  selling  point  for  MHD  propul¬ 
sion  systems  for  Naval  propulsion  application.  Therefore,  the  principal  starting  point  for 
this  project  in  fiscal  year  (FY)  1990  will  be  to  study  theoretically  the  physical  principles 
involved  in  the  generation  of  acoustic  noise  in  MHD  propulsors  on  submerged  ve¬ 
hicles.  Little  reliable  scientific  information  exists  on  acoustic  signature  of  these  MHD 
systems. 

GOALS  OF  THIS  TECHNICAL  REPORT 

The  MHD  propulsion  program  at  DTRC  was  initiated  late  in  FY89  with  a  small 
theoretical  effort  to  investigate  MHD  propulsion  efficiencies,  and  to  develop  plans  for 
possible  MHD  experiments  to  study  acoustic  noise  in  later  work.  This  technical  report 
documents  the  initial  theoretical  effort  on  the  MHD  program. 

A  literature  search  was  performed  on  marine  MHD  propulsion.*-17  Since  the  nu¬ 
merical  results  in  the  literature  were  often  too  incomplete  and  sketchy  for  reliable 
decision-making,  a  fundamental  set  of  equations  which  describe  the  propulsive  efficien¬ 
cies  of  an  ideal  rectangular  duct  propulsion  system  was  developed  into  a  computer 
program.  Assumptions  in  the  mathematical  model  were  generally  made  which  tend  to 
maximize  the  propulsive  efficiencies  of  the  MHD  system.  These  numerical  results  can  be 
used  for  engineering  estimates  of  the  propulsive  efficiencies  of  real  MHD  propulsion  sys¬ 
tems.  The  computer  program  was  exercised  for  a  wide  variety  of  cases;  the  work  is 
described  in  this  report. 

This  report  begins  with  a  background  description  of  fundamental  propeller 
theory.  It  is  shown  that  the  propeller  propulsion  theory  is  approximately  analogous  to  a 
MHD  propulsive  rectangular  duct.  Next,  the  propulsion  theory  of  an  MHD  rectangular 
duct  was  derived  from  the  basic  principles  of  fluid  dynamics.  From  this  work,  a  oue- 
dimensional  control  volume  mathematical  model  for  an  ideal  rectangular  MHD  propul¬ 
sive  duct  was  developed,18’ 19  and  the  numerical  results  for  propulsive  efficiency  of 
different  system  conditions  are  presented  here  in  the  form  of  curves  and  tables. 

The  MHD  propulsion  mathematical  model  will  be  modified  in  FY  90  to  include 
channel  end  effects.  The  theory  of  Hughes  and  McNab20  was  developed  for  MHD  liquid 
metal  pumps  and  will  be  applied  to  MHD  seawater  propulsion  channels.  It  is  expected 
that  these  end  effects  will  degrade  the  propulsive  efficiencies  calculated  for  the  ideal 
MHD  propulsion  model  without  end  effects.  End  effects  are  known  to  have  a  degrading 
effect  on  MHD  liquid  metal  pumps.  However,  channel  duct  design  in  MHD  propulsion 
systems  can  be  used  to  lessen  these  end  effects. 
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TRADITIONAL  PROPELLER  THEORY  AND  APPROXIMATE 
EXTENSION  TO  MHD  DUCTS 


BACKGROUND 

We  will  begin  the  discussion  with  the  traditional  propeller  theory,  since  this  topic 
closely  relates  to  the  understanding  of  the  magnetohydrodynamic  (MHD)  propulsion 
problem.  The  ship  being  considered  is  driven  at  a  constant  velocity  V  with  a  single  pro¬ 
peller.  The  frame  of  reference  is  chosen  to  be  moving  with  the  ship.  Thus,  the  fluid 
unaffected  by  the  propeller  has  velocity  V  and  pressure  PD.  The  propeller21  is  simply 
modelled  as  a  porous  disk  with  area  Ap  through  which  water  flows.  The  water  flowing 
through  the  disk  experiences  a  pressure  increase,  A P,  as  it  passes  through  the  disk  repre¬ 
senting  the  propeller.  We  will  designate  the  pressure  on  the  upstream  side  of  the  disk  as 
Pu,  which  is  less  than  the  ambient  pressure  Pa.  Because  the  relationship  Pu  <  P0  holds, 
the  water  is  accelerated  from  V  to  Vp  as  it  approaches  the  propeller.  Because  of  this  fluid 
acceleration,  the  water  passing  through  the  propeller  comes  from  a  far  upstream  stream 
tube  with  area.  A,  which  is  larger  than  Ap.  As  the  water  passes  through  the  disk,  its  ve¬ 
locity  remains  Vp,  by  conservation  of  mass,  but  its  pressure  increases  to  (Pu  +  AP),  which 
is  greater  than  the  ambient  pressure  P0.  The  high  pressure  water  that  leaves  the  propeller 
must  merge  with  the  surrounding  water  at  P0,  so  the  water  accelerates  to  the  exit  velocity 
Ve  which  is  greater  than  the  velocity  Vp.  There  is  a  corresponding  reduction  in  exit  area 
to  At  and  a  return  of  the  pressure  to  P0.  In  the  moving  reference  frame  far  downstream  of 
the  propeller,  we  have  a  jet  with  elevated  velocity'  V,  and  cross-sectional  area  A, ,  sur¬ 
rounded  by  water  with  velocity  V  and  pressure  at  PQ.  Figure  1  shows  a  diagram 
representing  this  traditional  propeller  system. 

There  are  four  stations  on  the  diagram.  Far  upstream  at  station  1 ,  the  ambient  pres¬ 
sure  is  Pc,  the  area  of  the  free  stream  tube  is  A,  and  the  velocity  of  the  fluid  is  V.  At 
station  2,  immediately  upstream  of  the  propeller,  the  pressure  Pu  is  less  than  the  ambient 
pressure  P0,PU<  P0,  the  area  of  the  propeller  is  Ap  and  the  velocity  into  the  propeller  is 
Vp.  At  station  3,  immediately  downstream  of  the  propeller,  the  pressure  jumps  from  Pu 
immediately  upstream  to  Pu  +A P  immediately  downstream.  Pu  +AP  is  greater  than  the 
ambient  pressure  (Pu  +AP)  >  P0,  the  area  of  the  propeller  is  Ap,  and  the  velocity  coming 
out  of  the  propeller  i sVD.  Far  downstream  the  ambient  pressure  is  P0  ;  the  exit  area  is  A,, 
and  the  exit  velocity  is  Vt. 

From  conservation  of  mass  the  following  relationship  holds:  VA  =  VDAp  =  VtAt . 
The  law  of  conservation  of  energy  can  be  used  to  represent  the  conservation  of  mechani¬ 
cal  energy  per  unit  of  mass.  This  water  flowing  from  a  to  b  along  a  streamline  can  be 
represented  by: 


Pa  .  V/  _  Pb  .  Vb 


—  +  — —  +  gia  =  —  +  — —  +  gzb  +  Ws  +  W/=  Constant  along  streamline  (1 ) 

Q  2  Q  2 


where: 


V2 

2 


vertical  coordinate,  with  positive  direction  upward, 
kinetic  energy  per  unit  mass  , 
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AXISYMMETRIC  FREE  STREAM  TUBE 


o 


Fig- 1*  Stream  tube  for  propeller  system. 


Ws  =  work  done  by  fluid  per  unit  mass, 

=  energy  lost  to  friction  flowing  from  a  to  b,  per  unit  mass, 
gz  =  potential  energy  in  elevation  per  unit  mass,  and 

—  =  potential  energy  in  pressure  per  unit  mass. 

e 

Like  all  potential  energies,  both  of  these  have  arbitrary  datum,  i.e.,  we  can  measure  all 
elevations  from  any  horizontal  plane  and  all  pressures  from  any  arbitrary  constant  pres¬ 
sure. 

Some  authors  refer  to  Eq.  1  as  a  modified  Bernoulli  equation.  The  following  as¬ 
sumptions  are  made  in  order  to  simplify  Eq.  1: 

•  Flows  are  essentially  horizontal,  and  we  can  subtract  out  hydrostatic  pres¬ 
sure  variation,  so  we  ignore  the  gz  terms. 


•  It  is  assumed  that  the  propeller  is  frictionless  (not  for  MHD  duct)  so  that 
Wf  =  0  for  the  propeller  case. 

Referring  to  Figure  1  for  the  diagram  for  the  propeller  flow,  Eq.  1  from  stations  1  to 


Po  ,  V2  =  Pu  |  Vp2 
Q  +  2  e  +  2 


(2) 
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Note  that  only  the  propeller  does  work  on  fluid  (negative  W*).  Therefore,  Ws  =  0  from 
stations  1  to  2,  and  from  stations  3  to  4: 


PU  +  AP  |  Vp2  =  Pe  ;  Ve2 
Q  2  Q  +  2 


Equation  1  from  stations  2  to  3  is: 

A P 

Ws  = - .  (4) 

Q 

We  will  define  T  here  to  represent  the  forward  thrust  on  the  propeller  disk  which  is 
equal  to  the  afterward  negative  thrust  of  propeller  on  water.  The  thrust  from  the  propeller 
is: 


T=(AP)Ad  .  (5) 

The  thrust  -  T  on  the  submerged  vehicle  is  caused  by  the  reaction  force  and  is  equal  in 
magnitude  and  opposite  in  direction  to  the  force  on  the  fluid  T. 

The  afterward  thrust  on  the  water  from  the  propeller  can  be  calculated  from  th-  con¬ 
servation  of  linear  momentum  for  steady,  incompressible  flow  through  a  control  volume 
(that  is,  a  fixed  volume  in  space  with  control  surface).  (See  Fig.  2.) 

T=YF  =  Q  f \w„dA, 

contrr1  (6) 

surface 


dA 


OUTWARD  UNIT 
NORMAL  TO 
CONTROL  SURFACE 

— ^  — > 

Vn  =  Vn 


6 


DTRC— 90/009 


where: 

Vn  =  outward  normal  component  of  velocity  at  a  point  on  control  surface  , 

dA  =  differential  element  of  control  surface, 

dQ  =  V„  dA  =  volumetric  flow  across  dA  per  unit  time,  and 

F  =  sum  of  all  external  forces  acting  on  water  which  is  instantaneously  inside 

control  surface,  including  pressure  forces  from  water  outside  control  surface. 

Figure  3  shows  the  control  volume  for  the  propeller  with  the  force  T  acting  on  the 
water. 

For  the  control  surface  for  the  propeller  we  assume  that  P  =  P0  over  entire  sur¬ 
face.  Therefore,  for  this  problem  pressures  give  zero  net  force.  This  is  the  traditional 
assumption,  which  is  not  strictly  true.  With  this  assumption  the  only  external  force  for 
the  system  is  the  propeller  thrust  T.  At  station  1,  V„  =  -  V  (fluid  inflow)  and  Vy  =  +  V, 

and  the  entrance  area  equals  A.  At  station  4,  V„  =  +  Ve  (outward  flow  of  fluid), 

Vy  -  +  Ve ,  and  exit  area  equals  A,.  At  the  stream  tube  surface,  V„  =  0.  Substituting 
these  results  into  Eq.  6  results  in  the  well  known  expression  for  the  propeller  thrust: 

T  =  g(V)(-V)A  +  6(Ve)(Ve)Ae  ,  (7) 


o 


Fig.  3.  Control  volume  lor  propeller  mathematical  model. 
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But  from  the  incompressibility  of  the  water  the  following  relationships  holds: 


Q  =  VA  =  VDAD=VeAe  ,  (8) 

where  Q  =  volumetric  flow  rate  (m3/s) .  Thus  Eq.  7  can  be  conveniently  expressed  as: 

T  =  eQ(Ve-V)  ,  (9) 

where  Q  =  mass  flow  rate  (mass  per  unit  time),  and  gQ(V,  -V)  =  change  in  linear  mo¬ 
mentum  per  unit  time . 

First,  we  consider  a  traditional  manipulation  of  Eq.  1-9.  Then  we  separate  indepen¬ 
dent  from  dependent  variables  in  these  equations.  Equating  Eqs.  5  and  9,  we  obtain: 

(A P)AD  =  gQ(Ve-V)  ,  (10) 


Since  Q  =  VdAd ,  this  gives 

A P 

- =  VD(Ve-V)  .  (11) 

Q 

Subtracting  Eq.  2  from  Eq.  3: 

(A P)  1  -  ,  1 

— -±  =  -(Ve2-V2)=-(ye  +  V)(V'-V)  .  (12) 

Q  2  2 

Equating  Eqs.  11  and  12,  we  obtain 

VD  =  ^(Ve  +  V)  .  (13) 

This  is  a  full  solution  of  the  energy  equation  and  the  momentum  equation. 

Now  consider  dependent  and  independent  variables  involved  in  the  problem.  Take 
ship  velocity  V  and  propeller  disk  area  as  the  only  two  independent  variables 

Ad  =  ( n/4)D 2  .  The  force  T,  required  to  move  the  ship  at  constant  velocity  V,  is  uniquely 
determined  by  V  and  the  hull  shape.  Consider  a  streamlined  submerged  vessel  with  only 
skin  friction  drag.  Then  the  following  equation  holds: 

Skin  friction  drag  =  T  =  ^qV^CjSw  ,  (14) 


where: 

Sw  =  wetted  surface  area  of  vessel; 

Cf  =  skin  friction  coefficient,  which  is  a  function  of  Reynold’s  number;  and 
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qVLv 

Re*.  = -  ,  where  U  equals  vessel  length,  g  equals  density  of  water,  and  ft  equals  vis- 

M 

cosity  of  water. 

The  Reynold’s  number  can  be  used  to  determine  the  skin  friction  coefficient  Cf .  If 

1  33 

the  flow  is  laminar  (Ret  <5  x  105)  then  Cf  can  be  represented  by  Cf  =  — '  .  At 

the  transition  point  where  the  Reynold’s  number  has  the  range  5  x  105  <  Ret  <  107 ,  Cf 

0.074  1740  , 

can  be  represented  by  Cf  =  — — -  — —  .  In  the  turbulent  region  (Ret  >  10  ) ,  the 

frictional  coefficient  Cf  can  be  represented  by  the  Prandtl-Schlichting  equation, 

0.455  1700 

/_  (log  10  Ret)2-58  ~  Ret 

Therefore,  once  we  specify  V,  the  required  T  is  known  from  the  skin  friction  drag, 

Eq.  14.  Since  At>  must  be  specified,  the  following  relationship  holds  from  Eq.  5: 


A P  = 


(15.1) 


(15.2) 


Using  Eqs.  12  and  15.2,  we  obtain: 


T  ■ 


From  Eq.  13,  VD  is  easily  obtained: 


Using  the  relationship  VA  =  VpAo  =  V,A,  from  the  incompressibility  of  the  fluid,  we  ob¬ 
tain: 


A  =  i+  i+SJs 

2  \  An 


y/2- 


(18.1) 
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Cf  Sw 

1+“1 — 
Ad 


-1/2 


(18.2) 


Everything  is  now  determined  from  V  and  Ap  plus  Cf  and  Sw- 

APPROXIMATE  RELATIONSHIP  OF  PROPELLER  TO  MHD  DUCT 

We  will  consider  a  simple  MHD  rectangular  duct  like  that  of  Fig.4.  Converging  or 
diverging  nozzles  attached  to  a  simple  duct  in  general  tend  to  degrade  the  propulsive  effi¬ 
ciency  of  the  duct  due  to  surface  friction  effects.  Therefore,  we  will  discuss  a  simple 
rectangular  duct  in  this  section,  as  opposed  to  a  more  complicated  configuration,  because 
at  this  point  we  are  interested  in  the  upper  limits  of  propulsive  efficiency. 

The  duct  problem  is  approximately  analogous  to  the  traditional  propeller  problem 

where  the  force  T  is  produced  by  the  J  x  B  Lorentz  body  force  rather  than  the  propeller 
force.  For  a  detailed  discussion  of  this  analogy,  the  reader  should  consult  Appendix  A. 


x 


FREE  STREAM 


TUBE 


V0 


Fig.  4.  Simple  MHD  rectangular  duct. 
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Equations  16, 17, 18.1,  and  18.2  give  the  expressions  for  Ve,  Vp,A  and  Ac  The  duct  di¬ 
mensions  are  length  L  =  (Ay) ,  distance  between  electrodes  w  =  (Ax) ,  distance  between 
insulators  h  =  (At),  and  cross  sectional  area  Ap  =  wh.  The  Cartesian  coordinate  system 
(x,y,z)  is  shown  in  Fig.  4,  where  x,y ,z  are  the  unit  vectors.  Throughout  the  duct,  we  ig¬ 
nore  end  effects.  It  is  well  known  that  duct  end  effects  degrade  propulsive 
efficiency.  The  velocity  of  fluid  flow  in  the  duct  is  a  constant  V  =  Vpy,  because  the  duct 

has  a  constant  cross-sectional  area.  The  electric  current  density  is  designated  as  j  =  Jx, 
the  electric  field  as  E  =  Ex  =  (<p / w)x ,  and  the  voltage  difference  between  electrodes  as 
<p  .  The  magnetic  field  is  represented  by  B  -  -  Bz,  the  Lorentz  body  force  as 
jxB  =  JBy,  and  the  “back”  electromotive  force  is  VxB  =  -  VpBx. 

The  conservation  of  energy  equation  from  the  entrance  to  the  exit  of  the  duct  is 

V entrance  ~  V exit,  ~ entrance  ~  ' exit 


A P 

- =  -Ws-Wf  , 

Q 


where: 


-Ws=; - =  work  done  on  a  unit  mass  by  electromagnetic  body  force,  JxB ,  as 

Q 


mass  travels  from  entrance  to  exit,  and 


=  frictional  losses  per  unit  mass. 


where: 


Dh  -  hydraulic  diameter  = 


4 Ap  2wh 


n  —  utuiuvuiv  wmiiivivt  —  —  , 

Per  ( w  +  h ) 

Per  =  perimeter  =  2(w  +  h) ,  and 

/  =  friction  factor,  which  is  a  function  of  Rep,  —  and  M. 

Dh 

oVpDh 

Rtp  =  — — —  =  Reynold’s  number  based  on  mean  duct  velocity  and  hydraulic 
diameter,  and 

e  =  characteristic  wall  roughness  size,  which  is  considered  a  property  of  the  wall 

material.  For  example,  for  cast  iron  the  value  of  e  =  0.00085  ft,  while  for  wrought 
iron,  its  value  is  0.00015  ft. 

—  =  dimensionless  roughness  size. 

Dh 
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M  =  BDh 


1/2 


=  Hartmann  number  based  on  hydraulic  diameter,  and 


o  =  electrical  conductivity  of  fluid. 


Combining  Eqs.  15  and  19,  we  obtain  the  electric  current  density  required  to  over¬ 
come  the  frictional  resistance  in  the  duct  and  provide  a  pressure  rise,  AP,  needed  to  drive 
the  ship  at  a  steady-state  condition. 


J 


2 B 


C/Sv 

Ad 


(20) 


skin  friction  drag  of  hull  +  internal  friction  of  duct  flow 


The  only  role  of  ship  length  is  in  the  relationship  of  Rei  to  Cf  and  in  Sw.  Now  we 
must  determine  what  voltage  and  power  is  needed  to  produce  this  J. 

From  the  ^-component  of  Ohm’s  law  22: 


J -a 


(21.1) 


or 


(21.2) 


Equation  21  gives  the  voltage  difference  between  plates  required  for  this  thrust.  The  pro¬ 
pulsive  efficiency  rj  is  represented  by  the  relationship: 


V 


TV 

tf>I 


(22) 


T  is  given  by  Eq.  5  and  <t>  is  given  by  Eq.  21 .2,  in  terms  of  J  and  VD.  The  total  cur¬ 
rent  across  the  rectangular  channel  is  /  =  JhL. 


12 
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MATHEMATICAL  MODEL  OF  MAGNETOHYDRODYNAMIC  (MHD) 
SHIP  PROPULSION  WITH  PROPULSIVE  EFFICIENCY 


FUNDAMENTAL  EQUATIONS  OF  MHD  PROPULSION 

This  section  presents  the  development  of  the  fundamental  equations  for  an  MHD 
propulsive  rectangular  duct  on  a  submerged  vehicle.  The  duct  configuration  under  con¬ 
sideration  is  presented  in  Fig.  5.  The  fluid  flow  into,  through,  and  out  of  the  duct  is 
represented  by  a  stream  tube  of  total  surface  area  2  =  2]  +  E2+23+A+A<  and  total 
volume  V  '=  V\  +  V2'  +  V3' .  The  duct  together  with  the  stream  tube  is  surrounded  by  an 
ambient  pressure  P0.  The  duct  is  of  uniform  cross-sectional  area  Ap.  The  entrance  area 
and  exit  area  of  the  stream  tube  are  A  and  Ae  respectively.  The  pressure  at  the  entrance  of 
the  stream  tube  is  P0 ,  at  the  duct  entrance  Pu ,  at  the  duct  exit  Pu  +  AP ,  and  at  the  stream 
tube  exit  P0 .  The  following  inequalities  hold  for  the  pressures:  P0  >  PM ,  Pu  <  Pu  +  AP , 
and  Pu  +  AP  >  P0.  The  velocity  of  the  fluid  entering  the  stream  tube  is  V,  entering  the 
duct  Vp,  leaving  the  duct  Vp,  and  leaving  the  stream  tube  Vt.  The  unit  normal  to  the 
stream  tube  n  is  outward  from  the  surface. 

The  Cartesian  coordinate  system  (x,y,  2)  has  the  origin  fixed  in  the  center  of  the 
duct  and  moves  with  the  duct  at  constant  velocity  V.  The  observer  is  therefore  considered 

to  be  travelling  with  the  duct.  The  constant,  homogeneous  magnetic  induction  B  is  in  the 
negative  z-direction.  The  constant,  homogeneous  current  density  vector  lines  are  in  the 
positive  x-direction.  The  fluid  flows  in  the  positive  y-direction. 


TUBE 


2  *  a  +  Aq  +  Si  +  y.p  +  2,3 
V  m  v/  +  Vz  +  V3 


Fig.  5.  MHD  rectangular  duct. 
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Using  the  above  rectangular  duct  configuration,  we  shall  develop  the  general  funda¬ 
mental  equations  for  deriving  the  propulsive  efficiency  of  the  rectangular  MHD  duct  in  a 
steady  state  condition.  The  following  fundamental  principles  from  theoretical  physics  for 
fluid  flow  will  be  used  for  the  derivation. 

•  Conversation  of  linear  momentum  for  a  stationary  control  volume 

•  Rate  of  change  of  kinetic  energy  for  a  stationary  control  volume  (Energy 
Theorem) 

•  The  principle  of  action  and  reaction 

•  Ohm’s  law  for  electrical  conductors 

This  derivation  implies  the  conservation  of  mass,  linear  momentum,  energy,  and  angular 
momentum  for  the  system. 

Also,  the  approximations  to  these  equations  are  developed  for  the  mathematical 
model  that  was  used  to  produce  the  numerical  computer  results  presented  in  the  report. 
This  development  follows  well-known  concepts  of  fluid  mechanics  for  fluid  flow  through 
a  stationary  control  volume  and  is  discussed  in  detail  to  clarify  the  concepts  involved  in 
an  MHD  duct. 

Conservation  of  Linear  Momentum 

Newton’s  Second  Law  states  that  the  sum  of  all  external  forces  acting  on  an  assem¬ 
blage  of  masses  equals  the  rate  of  change  of  the  total  linear  momentum  of  the 
assemblage.  Our  assemblage  is  the  mass  of  fluid  instantaneously  inside  the  control  vol¬ 
ume  V '  in  Fig.  5.  Since  fluid  is  continuously  entering  and  leaving  our  fixed  control 
volume,  there  are  two  different  assemblages  of  masses  at  two  different  times.  A  change  of 
the  linear  momentum  of  the  assemblage  instantaneously  inside  the  control  volume  ap¬ 
pears  in  two  forms.  First,  the  velocity  at  points  inside  the  control  volume  may  be 
increasing  or  decreasing  with  time,  so  that  the  assemblage  of  masses  inside  the  control 
volume  at  time  (t+dt)  has  more  or  less  linear  momentum  than  the  assemblage  at  r.  These 
two  assemblages  differ  only  by  the  control  surface,  Z  =  A  +  At  +  Z  i  +  2  2+  2  3,  during 
the  short  interval  from  t  to  (t+dt).  This  first  change  of  linear  momentum  is  given  by 

-ff  QV'dV' ,  (23.1) 

ot  Jv  ■ 

where  V;  are  the  components  of  the  velocity  V,  and  dV’  is  a  differential  volume 
element.  The  second  realization  of  a  linear  momentum  change  is  that  the  fluid  leaving  the 
control  volume  may  have  more  or  less  linear  momentum  that  the  fluid  entering  ii  This 
second  change  of  linear  momentum  is  given  by 

j^QViVknkdS,  (23.2) 
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where  nk  are  the  components  of  the  outward  unit  normal  n  at  each  differential  element  dS 
of  the  control  surface  5.  The  volume  of  fluid  crossing  the  element  dS  of  the  control  sur¬ 
face  per  unit  time  is  Vk  nk  dS,  which  is  positive  or  negative  for  fluid  leaving  or  entering, 
respectively,  while  the  linear  momentum  of  this  fluid  per  unit  volume  is  gV, .  Therefore, 
Newton’s  Second  Law  for  a  Control  Volume,  which  is  an  Open  System,  is 

2^  =  t[  QVjdV'+f  QV,vknkdS  ,  (24) 

dt  J  v  •  Jz 

where  is  the  sum  of  all  external  forces  acting  on  the  mass  instantaneously  inside 
the  control  volume.  The  sum  of  forces  is  given  by 

%Fi=  (bidV’+f  ijdS  ,  (25.1) 

where  b,  is  the  external  body  force  per  unit  volume  at  each  point  inside  the  control  vol¬ 
ume,  and  r,  is  the  external  traction  force  per  unit  area  at  each  point  on  the  control  surface. 
For  a  viscous  fluid 

ti  =  -Pni  +  Ofc  nk  ,  (25.2) 

where  P  is  the  pressure  and  a*  is  the  viscous  shear  stress.23  If  P  represents  the  deviation 
of  the  actual  pressure  from  the  hydrostatic  pressure,  than  we  can  ignore  the  gravitational 
body  force,  so  that  our  only  body  force  is  the  Lorentz  force  per  unit  volume, 

bi  =  €ijk  JjBk  ,  (25.3) 

where  ei;*  is  the  permutation  tensor,  Jj  is  the  electric  current  density,  and  Bk  is  the  mag¬ 
netic  induction.  We  consider  a  steady  flow,  so  that  the  velocity  at  each  point  inside  the 
control  volume  is  independent  of  time,  and  Newton’s  Second  Law  is 

6ViVknkdS  =  eijkJjBk  dV’+  dS-  j  Pm  dS  .  (26) 

The  sum  of  the  external  forces  on  the  right-hand  side  act  on  the  seawater  which  is 
inside  the  control  surface  S.  The  first  term  is  the  total  Lorentz  body  force  on  the  fluid  in 
the  rectangular  channel  plus  end  effects.  The  second  term  represents  the  frictional  forces 
of  the  fluid  along  the  surfaces  2  x  +  2 ,  +  2  3.  The  frictional  force  along  the  inside  walls 

of  the  duct  2  2  is  more  important  than  the  frictional  force  along  2  j  +  2  3.  The  third 

term  includes  the  pressure  forces  from  the  fluid  external  to  the  control  surface.  The  fol¬ 
lowing  integrals  can  be  evaluated  in  Eq.  26. 
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j  i ?  ViVk  nkdS  =  gQ(Ve -V)x  (27.1) 

where  Q  =  VA-  VpAp  =  VfAt. 


I 


PndS  =  0 


(27.2) 


j"  Oik'nk  dS  =  j  Oiknk  dS  (27.3) 


Substituting  these  integrals  in  Eq.  26  results  in  the  following  equation: 


eQ(.ve-v)i  = 


(7x  B)dV'  + 
v 


dS  . 


For  the  interior  surface  of  the  duct  2  ,, 


(28.1) 


Oik  rijfc  =  T„, 


~^feVD2x  , 


(28.2) 


where  rw  is  the  viscous  wall  shear  stress  and/is  the  Moody  friction  factor.  If  we  ignore 
all  electromagnetic  end  effects,  the  ^-component  of  Eq.  28.1  is 


QQ(Ve-V)  =  JBV2'-^fQ  VD2L  (Per) 

T  =  Lorentz  force  -  Internal  Friction  force 
V2'  =  LWh 


where  L,  W,  h,  VV>  and  (Per)  =  2(w+h)  are  the  length,  width,  height,  interior  volume, 
and  cross  sectional  perimeter,  respectively.  Equation  28.3  is  also  obtained  by  substituting 
Eq.  19  into  Eq.  10.  The  term  T  is  designated  as  the  effective  thrust  on  the  fluid  down  the 
channel.  Equation  28.3  is  the  fundamental  momentum  equation  for  the  operation  of  the 
MHD  duct. 

Principle  of  Action  and  Reaction 

The  reaction  force-T  to  the  effective  Lorentz  body  force  T  propels  the  submerged 

vehicle  at  a  steady  state  velocity  V  in  the  fluid.  The  reaction  force-T  balances  the  hydro- 
dynamic  resistance  R  of  the  vehicle.  The  reaction  force-7  is  caused  by  the  magnetic 
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induction  B  from  the  current  density  in  the  channel  interacting  with  the  current  in  the 
channel  magnetic  coil  electrical  system.  Therefore,  the  useful  power  in  propelling  the 
submerged  vehicle  by  the  MHD  duct  is  expressed  as 

Useful  power  in  propelling  vehicle  =  TV  (29) 

The  reaction  to  the  force  JBVf  acts  directly  on  the  MHD  channel  magnets. 

For  a  more  detailed  discussion  of  the  reaction  force  see  Appendix  A. 

Rate  of  Change  of  Kinetic  Energy  (Energy  Theorem) 

The  dot  product  between  Newton’s  Second  Law  for  a  single  mass  and  the  mass’s 
velocity  vector  gives  an  equality  between  the  rate  of  change  of  the  mass’s  kinetic  energy 
and  the  work  done  on  it  per  unit  time.  When  we  sum  the  Newton’s  Second  Law  for  an 
assemblage  of  masses,  all  internal  forces  between  the  masses  sum  to  zero,  since  the  forces 
between  two  masses  are  equal  and  opposite.  When  we  sum  the  kinetic  energy  equation 
for  an  assemblage  of  masses,  the  work  done  by  the  internal  forces  does  not  sum  to  zero. 

In  a  continuum  model,  the  internal  forces  are  represented  by  the  pressure  and  viscous 
stresses.  The  pressure  force  is  conservative  for  a  liquid,  so  it  can  only  transfer  kinetic  en¬ 
ergy  between  masses  with  no  net  change  in  the  total  kinetic  energy  of  the  assemblage. 

The  viscous  stresses  are  not  conservative  and  produce  an  irreversible  conversion  of  kinet¬ 
ic  energy  to  heat.  The  rate  of  energy  conversion  per  unit  volume  is  called  the  viscous 
dissipation  0  and  is  proportional  to  the  product  of  the  viscous  stress  between  neighbor¬ 
ing  masses  and  their  difference  in  velocity,  as  reflected  in  the  velocity  gradient, 

dV 

0  -Oik — -  =  viscous  dissipation.  (30) 

dxk 

The  rate  of  change  of  the  kinetic  energy  of  the  assemblage  of  masses  instantaneously  in¬ 
side  the  control  volume  is  equal  to  the  rate  of  increase  of  kinetic  energy  inside  the  control 
volume,  plus  the  difference  between  the  kinetic  energies  in  the  fluid  entering  and  leaving 
the  control  volume.  Thus,  the  rate  of  change  of  kinetic  energy  for  a  control  volume  is 

4 ’  f  d/2 )QV2dV'+  f  (1/2 )gV2Vk  nk  dS 

dt  Jv  JZ 

=  f  (biVi)  dV'+  f  tiVi  dS-  f  Ouf^dV'  ,  (31) 

Jv  Jz  Jv  dxk 

where  again  b,  and  r,  represent  the  external  body  force  and  surface  traction. 

Alternatively,  we  can  begin  with  the  First  Law  of  Thermodynamics,  which  states 
that  the  heat  added  to  an  assemblage  of  masses,  plus  the  work  done  by  the  external  forces, 
equals  the  change  in  total  energy,  which  is  the  sum  of  the  kinetic  and  internal  energies. 
The  internal  energy  represents  the  internal  energy  of  the  molecular  states.  We  can  split  the 
First  Law  of  Thermodynamics  into  two  separate  equations,  with  the  heat  and  work  in  the 
internal  and  kinetic  energy  equations,  respectively.  In  making  this  split,  we  must  add 
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equal  and  opposite  source  terms  to  each  equation  to  represent  the  energy  transfer  between 
internal  and  kinetic  energies.  In  a  gas,  internal  energy  can  be  converted  to  kinetic  energy 
through  volumetric  expansion,  but  in  a  liquid,  the  only  energy  transfer  is  the  viscous  con¬ 
version  of  kinetic  energy  to  internal  energy.  After  the  split  of  the  First  Law,  the  kinetic 
energy  equation  is  Eq.  31. 

Introducing  Expressions  25.2  and  25.3,  where  again  P  represents  the  deviation  from 
hydrostatic  pressure,  Eq.  31  becomes 

—  f  (1/2 )gV2dV'+  f  [(1/2 )gV2  +  P)  Vknk  dS 
dt  Jv  Jl 

=  J  (vfiijtJjBk)  dV'+  j nk  dS  -  aik  dV'  .  (32) 

Our  flow  is  steady,  so  that  the  first  term  in  Eq.  32  is  always  zero.  We  can  apply 
Eq.  32  separately  to  the  smaller  control  volume  V/  or  vy  in  Fig.  5.  We  are  neglecting 

viscous  stresses  in  these  volumes,  and  JxB  =  0  here  since  we  neglect  electromagnetic 
end  effects.  Therefore,  Eq.  32  becomes 

Q[(1/2)qVd2  +  PM  -  (l/2)gV2  -P0]  =  0  (33.1) 

Q[(\/2)gV2  +  P0  -  (1/2 )qVd2  -P/i-AP]  =  0  (33.2) 

for  volumes  Vj'  and  V3',  respectively.  The  Eqs.  33.1  and  33.2  are  identical  to  the  con¬ 
servation  of  energy  equations  ( Eqs.  2  and  3  )  after  dividing  through  by  Q.  When  we 
apply  Eq.  32  to  the  other  smaller  control  volume  V2  in  Fig.  5,  we  must  determine  the 
two  viscous  terms  on  the  right-hand  side.  We  use  a  one-dimensional  flow  model  with 

_  ^  Wi 

V  =  Vpx  throughout  V2.  Therefore,  — -  s  0,  and  the  viscous  work  is  the  work  done 

dxk 

by  the  fluid  moving  at  velocity  Vpx  against  the  viscous  wall  shear  stress,  ~twx  at  X2, 
which  contributes 


-  xwVDUPer)  =  -  UQVDlL{Per)  (34) 

O 

to  the  right  side  of  Eq.  32.  In  reality,  the  fluid  at  2  2  must  be  at  rest,  so  that 

Oik'Vflk  -0  at  2 it  >s  the  final  viscous  dissipation  term  which  represents  the 

work  done  against  viscous  forces.  Actually,  this  distinction  is  moot  because  the  Moody 
friction  factor  is  defined  by  equating  the  Expression  34  to  the  power  input  needed  to 
maintain  a  constant  kinetic  energy  in  spite  of  frictional  losses.  For  the  control  volume  for 
VY,  Eq.  32  becomes 
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(35) 


Q(AP)  =  VDJBV2'~-fg V^UPer)  . 

O 

If  we  divide  this  equation  by  gQ  =  qV^Ad,  it  becomes  Eq.  19. 

Electromagnetic  Energy  Input 

In  our  frame  of  reference  moving  with  the  duct,  there  is  an  electric  filed  £,  a  mag¬ 
netic  induction  B ,  and  an  electric  current  density  7  in  the  fluid  with  electrical 
conductivity  o  inside  our  control  volume.  The  electromagnetic  energy  input  to  the  fluid 
per  unit  time  is24' 


If  we  introduce  Ohm's  law  without  Hall  effects, 


7  =  a  [e  +  V  x  £)  , 


in  order  to  eliminate  £,  then  Expression  36  becomes 


j  (J2/o  )dV'+  |  V-(j  X§)  dV'  . 


The  second  term  in  Eq.  38  is  the  work  done  on  the  fluid  by  the  Lorentz  body  force  (Eq. 
25.3),  and  this  term  appears  in  our  kinetic  energy  equation  (Eq.  32).  The  first  term  in  Eq. 
38  represents  a  heat  input  because  an  electric  current  flowing  in  a  conductor  produces 
heat.  This  heat  input  is  called  the  Joulean  heating  and  appears  as  a  term  in  the  conserva¬ 
tion  of  internal  energy  equation,  along  with  the  viscous  dissipation  (Eq.  30)  as  two 
volumetric  heat  sources.  Since  we  cannot  get  any  useful  work  out  of  the  internal  energy 
of  a  liquid,  the  Joulean  heating  represents  an  energy  loss  for  our  system.  Equation  36  rep¬ 
resents  the  total  electrical  power  input  to  the  MHD  system  and  is  equal  to 

JEV2'  =I<f>  =  VDJBV2'+I2Re  ,  (39) 


where  /  =  JLh  =  the  total  current  between  the  electrodes,  <p  =  Ew  =  total  voltage  differ¬ 
ence  between  electrodes,  h  =  height  between  insulating  walls,  w  =  width  between 
electrodes,  and  Rt  =  ( w/oLh )  is  the  electrical  resistance  of  the  duct. 
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Efficiencies  For  MHD  Propulsion 

The  basic  propulsion  equations  for  the  MHD  rectangular  duct  were  derived  and 
presented  previously  in  the  report.  We  assume  that  the  following  quantities  are  known: 

V  =  velocity  of  the  vessel; 

Cf  =  skin  friction  coefficient  of  the  vessel; 

S*  =  wetted  surface  area  of  the  vessel; 

L,  h’,  h  =  length,  width  and  height  of  MHD  propulsion  duct; 

B  =  strength  of  uniform  magnetic  induction  inside  duct. 

Other  quantities  are  given  by  the  following  equations: 

T  =  skin  friction  drag  resistance  of  vessel  which  equals  thrust  developed  by  propulsion 
unit,  from  Eq.  14; 

A P  =  pressure  rise  through  the  duct,  from  Eq.  15; 

Ve  =  exit  velocity  from  control  volume,  from  Eq.  16; 

Vo  =  velocity  inside  duct,  from  Eq.  17; 

A  =  area  at  the  entrance  of  the  control  volume,  from  Eq.  1 8.1 ; 

Ae  =  area  at  the  exit  of  the  control  volume,  from  Eq.  18.2; 

J  =  uniform  electric  current  density  between  electrodes  inside  duct,  from  Eq.  20; 

<p  =  total  voltage  difference  between  electrodes,  from  Eq.  21 .2; 
t]  =  propulsive  efficiency,  from  Eq.  22. 

For  the  purposes  of  interpretation,  we  can  split  the  propulsive  efficiency  into  three 
separate  efficiencies, 


V 


TV 

7* 


=  nsund  ■ 


(40.1) 


The  electrical  efficiency 


VD  JBLwh  I2Re 
Ve~  (pi  "1_  (pi 


(40.2) 


is  the  ratio  of  the  work  done  by  the  Lorentz  force  on  fluid  per  unit  time  to  the  total  electri¬ 
cal  power  input  to  duct.  Equation  39  has  been  substituted  in  order  to  get  the  second 
expression,  which  shows  that  (1  - r\t )  is  the  fraction  of  the  electrical  input  power  which  is 
lost  to  Joulean  heating.  The  thrust  efficiency 


TVp  {\/%)fQV[?L{Per) 

VD  JBLwh  VD  JBLwh 


(40.3) 


is  the  ratio  of  the  energy  gained  by  the  fluid  per  unit  time  inside  the  duct  to  the  work  done 
by  the  Lorentz  force  per  unit  time.  Equation  28.3  has  been  substituted  in  order  to  obtain 
the  second  expression,  which  shows  that  (1  -rjt)  is  the  fraction  of  the  work  done  by  the 
Lorentz  force  which  is  lost  to  the  internal  friction  of  the  duct.  The  hydraulic  efficiency 
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_  TV  d/2 )0Q(Ve-V)2 

Vd~  TVd  TVd 


(40.4) 


is  the  ratio  of  the  power  driving  the  vessel  to  the  energy  gained  by  the  fluid  inside  the 
duct  per  unit  time.  Equations  9  and  13  were  substituted  to  obtain  the  second  expression. 

In  a  fixed  frame  of  reference,  the  water  entering  our  control  volume  has  zero  velocity, 
while  that  leaving  has  a  velocity  (Ve-V).  The  kinetic  energy  in  the  stream  of  the  duct  rep¬ 
resents  a  loss  of  energy  per  unit  time,  which  is  (1  /2)gQ(Ve  -  V)2 . 

Up  to  this  point  we  have  implicitly  assumed  that  the  MHD  propulsion  duct  is  inside 
the  vessel’s  hull,  so  that  the  only  skin  friction  drag  resistance  is  that  of  the  hull.  If  the  duct 
is  housed  in  an  external  pod  which  is  connected  by  a  strut  to  the  hull,  then  the  thrust  de¬ 
livered  to  the  hull  is  reduced  by  the  skin  friction  drag  on  the  pod  and  strut.  Thus,  Eq.  28.3 
must  be  replaced  by 

T  =  JBLwh  -  (1/8 )fQVD2L{Per)  -  ( 1  /2)gV  2CfpSp  ,  (41 ) 


where  T  is  the  thrust  delivered  to  the  hull,  which  balances  the  skin  friction  drag  of  the 
hull,  as  given  by  Eq.  14,  Cfp  is  the  skin  friction  coefficient  for  the  pod  and  strut,  and  Sp  is 
the  wetted  surface  area  of  the  pod  and  strut.  Then  Eq.  20  is  replaced  by 


T 

hJL 

BLwh 

2 B 

■VDl  + 


2  ■  CfP$P 


As  a  conservative  estimate  of  the  pod’s  skin  friction  drag,  we  use 


Cfp  =  0.075  log  io 


/  gVL  \ 
\  P  ) 


2 1  +0.0005  , 


(43.1) 


Sp  =  L(Per)  , 


(43.2) 


where  Cfp  is  given  by  a  variation  of  the  Schoenherr  equation  which  applies  for 
0.0026  S  Cfp  S  0.0037.  Of  course,  Sp  must  be  much  larger  than  L(Per)  because  the 

pod  must  contain  the  large  magnets  needed  to  produce  B.  Therefore,  Eq.  42  becomes 


/  =  — —  +  — —  f fVD2  +  4CfnV : 
BLwh  2BDh  V  fp 
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In  addition,  Eq.  17  is  replaced  by 


_27_+4C££l\ 
qAdV2  Dh  ) 


1/2 


(45) 


Equations  21.2, 44,  and  45  now  give  (p ,  and  Eq.  22  gives  the  propulsive  efficiency. 
Rather  than  apply  this  analysis  to  specific  submerged  vessels,  we  will  replace  Eq.  14  by 

T  -  KV2  ,  (46) 


with  K  =  1000,  2000,  3000, 4000  or  5000  kg/m  for  some  typical  submerged  vessels.  If  the 
MHD  propulsion  duct  is  in  an  external  pod,  then  we  use  Eqs.  21 .2, 43.1 , 44, 45,  and  46  to 
compute  the  propulsive  efficiency  which  we  denote  by  172-  If  the  duct  is  inside  the  ves¬ 
sel's  hull,  we  use  the  same  equations  with  Cfp  =  0  and  denote  the  propulsive  efficiency  by 
Vi- 


Numerical  Results 

The  objective  at  this  point  of  the  work  was  to  calculate  the  propulsive  efficiencies 
rji  and  172  of  the  one-dimensional  control  volume  MHD  system  at  different  steady  state 
conditions.  These  numerical  curves  presented  herein  will  determine  the  approximate  up¬ 
per  limits  of  the  MHD  propulsive  vehicular  system  efficiencies,  which  can  be  used  for 
engineering  judgments  in  estimating  the  propulsive  efficiencies  under  different  conditions 
of  “real  world”  systems  used  on  submerged  vehicles.  The  efficiency  is  defined  as  the  use¬ 
ful  power  to  propel  the  vehicle  divided  by  the  total  electrical  power  to  the  duct  I<j) , 
where  I  is  the  total  current  across  the  duct  and#  is  the  voltage  drop  across  the  duct.  *  For 
the  efficiency  17 1  the  useful  power  to  propel  the  duct  is  defined  mathematically  as  the 
ideal  thrust  in  the  duct  multiplied  by  the  vehicular  speed  V  minus  the  frictional  force  in 
the  duct  (1  /8 )\fgVD2L(Per)\  multiplied  by  V.  The  efficiency  r]2  corresponds  to  an  MHD 
duct  in  the  form  of  a  pod  attached  by  struts  to  the  vehicle  hull.  The  propulsive  efficiency 
172  is  defined  analogously  to  the  propulsive  efficiency  rj\ .  In  this  case,  the  useful  power 
to  propel  the  submerged  vehicle  is  the  ideal  thrust  in  the  duct  multiplied  by  V  minus  both 
the  frictional  force  in  the  duct  due  to  fluid  flow  {l/S)\fgVp2L(Per)]  and  the  force  d-e  to 
the  fluid  flow  outside  the  pod  (1/2 )[gV2L(Per)Cfp] ,  where  both  frictional  terms  are  mul¬ 
tiplied  by  V.  From  these  definitions  rj\>rj2  under  the  same  propulsion  conditions. 

The  MHD  propulsive  mathematical  model  using  equations  from  the  previous  sec¬ 
tion  were  programmed  on  the  computer.  All  equations  are  in  the  SI  system  of  units.  The 
MHD  rectangular  duct  system  was  considered  to  be  propelling  a  submerged  vehicle  at  a 
steady-state  condition  V.  Transient  conditions  were  not  considered  at  this  point  of  the 

*  Note  that  this  number  does  not  account  for  power  generation  or  distribution  efficiencies. 
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work.  End  effects,  bubble  generation  power  losses,  and  the  effects  of  bubbles  in  the  duct 
were  also  omitted.  Computer  calculations  were  performed  using  various  input  parameters 
for  the  rectangular  duct  and  various  hydrodynamic  resistances  for  the  submerged  vehicle. 

To  properly  characterize  the  MHD  propulsive  system  propelling  a  submerged  ve¬ 
hicle  of  hydrodynamic  resistance  3,000  V2  newtons,  we  present  the  base  line  numerical 
numbers  in  Table  1.  We  consider  a  rectangular  duct  height  h  of  4  meters,  width  w  of 
4  meters,  and  length  L  of  10  meters  (see  Fig.  6).  A  realistic  magnetic  induction  B  of  6  tes¬ 
la  across  the  duct  was  chosen.  The  electrical  conductivity  o  and  density  of  seawater  g 
were  chosen  to  be  4  siemens/m  and  1026  kg/m3,  respectively.  The  remainder  of  the  calcu¬ 
lated  data  with  dimensional  units  is  shown  in  Table  1.  The  Moody  friction  factor /( i.e., 
/=  (4Cfj)  for  inside  the  MHD  duct  corresponds  approximately  to  that  for  smooth  pipes 
for  all  numerical  results  in  this  report. 

Figure  7  presents  the  total  power  consumption  of  the  vehicular  system  versus  ve¬ 
hicle  speed  V.  The  propulsive  efficiency  rj  i  was  used  in  this  plot.  The  solid  curve 


Fig.  6.  Propulsive  magnetohydrodynamic  rectangular  duct  configuration. 
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Table  1.  Baseline  numerical  results  for  square  rectangular  duct. 


Height  of  square  duct,  h 

4  meters 

Joulean  heat  loss  in  duct, 

I2  R 

4.59  x  1 06  watts 

Width  of  square  duct ,  w 

4  meters 

Electrical  current  in  duct,  / 

1.36  x  104  amperes 

Length  of  square  duct  L 

10  meters 

Coefficient  of  skin  friction 
outside  in  duct,  Ctd 

0.00262 

Magnetic  induction  of  duct,  B 

6  tesla 

Coefficient  of  skin  friction 
outside  of  duct,  Ct 

0.002637 

Velocity  of  vehicle,  V 

10  meters/second 

Total  input  power  to  propul¬ 
sive  system 

8.12x1 06  watts 

Velocity  of  fluid  in  duct,  VD 

10.8  meters/second 

Propulsive  efficiency,*?  1 

0.37 

Viscosity  of  seawater,  r\' 

1 03  newton-seconds/ 
meter 

Useful  power  to  propel  ve¬ 
hicle 

3.0  x  106  watts 

Electrical  conductivity  of 
seawater 

4.0  siemens/meter 

Power  to  overcome  skin  fric¬ 
tion  of  duct 

028  x  106  watts 

Density  of  seawater,  g 

1026  kilograms/meter3 

Power  to  overcome  external 
skin  friction 

0.0000 

Electrical  potential  across 
electrodes 

599.0  volts 

Power  lost  to  kinetic  energy 
aft 

025  x  106  watts 

Electrical  resistance  of 
duct,  R. 

0.025  ohms 

Hydrodynamic  resistance  of 
vehicle 

3.0  x  10s  newtons 

Electrical  current  density,  J 

338.9  amperes/meter2  1 

Nondimensional  parameters  of  duct  based  on  hydraulic  diameter,  Dh 
(Transverse  dimension  related  to  fully  developed  fluid  flow) 


force 

„  ijt  „  „  ^  _  ,-r  induced  magnetic  induction 

Magnetic  Reynold  s  number  Rem  =  VpDpjp  =  21.7  X  10~3  = - - - 

external  magnetic  induction 


Reynold’s  number 


„  oVnDh  .  M  inertial 

Rez>  =  ■=—7—  =  44.3  x  106  =  — - 


VISCOUS 


Interaction  number 


Hartmann  number 


Nd  _  oB2Ph  _  5  2  x  10 "2  =  Pondermotive  force 


QVd 


( _  v/2 

Md  -  LB\~\  =  3.8  x  10?  = 


inertial  force 

pondermotive  force 
viscous  force 


Nondimensional  parameters  of  duct  based  on  length,  L 
(Axial  dimension  related  to  boundary  growth  at  entrance  of  channel) 


Rez,  =  ^—7—  =  1 10.8  X  106 
V 

Remz.  =  VoLop  =  54.3  x  10"5 


oePl 

Nl  =  — — —  =  13  x  10-2 


QVd 


Ml 


5  x  103 
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presents  the  power  in  watts  required  to  propel  the  submerged  vehicle  at  constant  speed. 
The  dimensions  of  the  duct  are  the  same  as  in  Table  1 .  The  magnetic  induction  B  was 
6  tesla.  The  dotted  line  shows  the  total  input  power  required  to  propel  the  vehicle  system. 
The  percentages  above  the  dotted  line  show  the  propulsive  efficiency  rj\ .  It  should  be 
carefully  noted  that  a  representative  propulsive  efficiency  of  a  submerged  vehicle  driven 
by  a  single  screw  rotation  propeller  could  be  considered  to  be  about  85%.  These  curves 
show  clearly  that  the  propulsive  efficiency  of  this  MHD  system  decreases  significantly 
with  vehicle  speed  V,  and  that  the  propulsive  efficiency  is  much  less  than  that  of  a  propel 
ler. 

Figure  8  presents  the  propulsive  efficiency  curves  rj\  and  tj2  of  a  square  MHD 
propulsive  duct  with  a  6  tesla  magnetic  induction  B  across  the  duct,  propelling  a  sub¬ 
merged  vehicle  at  a  steady  state  speed  V.  The  duct  has  the  same  linear  dimensions  as  the 
duct  in  the  previous  figure.  The  hydrodynamic  resistance  of  the  vehicle  was  again  3,000 
V2  newtons  and  the  power  required  to  propel  the  vehicle  is  thus  3000  V3  watts.  The  pro¬ 
pulsive  efficiency  rj\>r}2,zs  would  be  expected  because  rj2  incorporates  the  loss  of 


DUCT  LENGTH,  L  -  10  m 

DUCT  HEIGHT,  h  -  4  m 

DUCT  WIDTH,  w  -  Am 

MAGNETIC  INDUCTION  IN  DUCT,  B  -  6T 

HYDRODYNAMIC  RESISTANCE  OF  VEHICLE,  R  -  3,000  V®  N 


VELOCITY,  V(m/s) 


Fig.  8.  Propulsive  efficiency  of  square  MHD  duct  versus  vehicle  speed  in 
meters  per  second. 
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power  due  to  the  flow  of  fluid  over  the  external  surface  of  the  duct.  This  corresponds  to 
the  duct  as  a  pod  attached  to  the  vehicle.  Both  curves  of  propulsive  efficiency  rj\  and  172 
show  a  significant  decrease  in  propulsive  efficiency  with  speed  V. 

Figure  9  presents  the  effective  thrust  T  in  newtons  versus  vehicle  speed  V  in  meters 
per  second  required  to  propel  a  submerged  vehicle  with  a  hydrodynamic  resistance  in  sea¬ 
water  of  order  3,000  V2  newtons.  The  thrust  curve  is  proportional  to  V 2  and  thus  rises 
steeply  with  V.  One  knot  is  equal  to  0.5144  m/second. 

Figure  10  presents  the  propulsive  efficiencies  rj  1  and  172  of  the  previously  de¬ 
scribed  propulsive  MHD  duct  versus  constant,  homogeneous  magnetic  induction  B  across 
the  duct  at  vehicular  speed  V  of  10  m/second.  Using  state-of-the-art  superconducting 
magnets,  a  magnetic  induction  of  6  to  10  tesla  in  a  large  volume  is  probably  the  highest 
magnetic  induction  that  could  be  obtained  in  practice.  The  6  and  10  tesla  limits  are  shown 

by  two  vertical  dotted  lines  on  the  figure.  The  propulsive  efficiencies  >71  and  tj2  increase 
significantly  with  the  magnitude  of  the  magnetic  induction  B.  The  power  to  propel  the 
vehicle  T  is  3000  V3  watts.  The  magnitude  of  the  separation  between  171  and  rj2, 

I?/!  -  J72I ,  increases  with  increasing  B. 

Figure  11  presents  the  propulsive  efficiencies  of  the  MHD  propulsive  duct  at  a  ve¬ 
hicular  speed  of  10  m/second  versus  electrical  fluid  conductivity  ff  in  a  6  tesla  magnetic 

induction  B.  The  propulsive  efficiencies  17 1  and  rj2  increase  significantly  with  fluid  elec¬ 
trical  conductivity  a  in  siemens  per  meter.  The  propulsive  efficiencies  of  the  MHD 
system  in  the  seawater  at  the  polar  ice  caps  (i.e.  o  =  0.4  siemens/m)  is  well  below  an  effi¬ 
ciency  of  0.10  for  either  171  or  rj2.  The  hydrodynamic  resistance  of  the  submerged 
vehicle  was  again  3000  V2  newtons. 

Figure  12  presents  the  propulsive  efficiencies  of  a  MHD  propulsive  duct  of  length  L 
equal  to  10  m  and  varying  square  cross-sectional  area  from  4  m2  to  approximately  81  m2. 
The  variation  of  cross-sectional  area  is  shown  in  the  horizontal  axis.  The  magnetic  induc¬ 
tion  across  the  duct  is  6  tesla,  and  the  power  required  to  propel  the  submerged  vehicle  is 

3000  V3  watts.  The  efficiencies  rj\  and  rj2  increase  significantly  with  increasing  cross- 

sectional  duct  area.  The  magnitude  of  separation  \rjx  -rj2\  between  rj\  and  172  becomes 
significant  at  large  duct  cross-sectional  areas.  It  should  be  realized  that  construction  of 
superconducting  magnets  for  a  duct  with  width  w  and  height  h  much  greater  than  4  m 
would  probably  be  impossible  with  present  technology. 

Figure  13  presents  the  propulsive  efficiencies  of  a  duct  of  width  w  of  4  m  and  a 
height  h  of  4  m  with  a  varying  length  L  in  meters  at  a  vehicular  speed  of  10  m/second. 

The  efficiencies  r)x  and  rj2  increase  slowly  with  duct  length  L  in  meters.  The  magnetic 
induction  across  the  duct  £  is  6  tesla.  The  maximum  duct  length  considered  in  this  figure 
is  20.0  m. 

Figures  14  through  18  show  the  effect  of  changing  the  MHD  rectangular  duct  width 
H’  through  the  values  1, 2, 4, 8,  and  16  m  on  the  duct  propulsive  efficiencies  rj\ .  On  each 
figure  for  a  specific  duct  of  width  tv  the  value  of  the  duct  height  h  is  increased  from  4  to 
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-  THRUST  REQUIRED  TO  PROPEL  SUBMERGED  VEHICLE 

AT  CONSTANT  SPEED  (THRUST  REQUIRED  TO  BALANCE 
HYDRODYNAMIC  RESISTANCE  OF  VEHICLE,  Ft) 

T  =  THRUST  REQUIRED  TO  PROPEL  VEHICLE 


3,000  V2  N 


VEHICLE  SPEED,  V(m/s) 


Fig.  9.  Thrust  required  to  propel  submerged  vehicle  at  constant  speed  V. 
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4.0  6.0  8.0  10.0  12.0  14.0  16.0  18.0  20.0 

MAGNETIC  INDUCTION  IN  DUCT  B  (T) 


Fig.  10.  Propulsive  efficiency  of  square  MHD  duct  versus  magnetic  in¬ 
duction  in  tesla. 

32  m.  The  MHD  duct  length  L  was  maintained  at  10  m.  On  each  figure  two  groups  of  nu¬ 
merical  curves  for  submerged  vehicle  speeds  of  2.5, 5,  and  10  m/second  are  presented  for 
magnetic  inductions  of  6  and  10  tesla.  The  vessel  hydrodynamic  resistance  was  3000  V2 
newtons.  The  corresponding  propulsive  efficiency  curves  rj\  on  each  figure  for  identical 
submerged  vehicle  speed  have  higher  tj  i  values  at  10  tesla  than  at  6  tesla  on  each  figure, 
as  would  be  expected.  On  each  figure  the  propulsive  efficiency  values  rj\  of  each  curve 
increases  with  duct  height  h  when  h  increases  from  4  to  roughly  10  m  and  then  decreases 
slightly  to  32  m.  The  main  point  that  each  figure  shows  is  that  the  corresponding  propul¬ 
sive  efficiency  curve  values  r)\  increase  with  duct  width  w. 

The  propulsive  efficiency  171  can  be  derived  from  Eq.  39  and  44  to  be 


TV 

=  — =  ■  BV  - 

I<f> 

h 


(l/S)fQVp2V2L  [l  +  (h/w) 


—  +  VdB 
Lha 
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ILSIVE  EFFICIENCY 


where  the  current  /  in  the  channel  is  determined  from 


Bw  4 B\  w  J 


In  Eq.  47,  the  frictional  term  containing  the  friction  factor/is  much  smaller  than 
BV,  and  in  Eq.  48  the  friction  term  containing /is  much  smaller  than  TIBw,  since  L  is  only 
10  m.  Increasing  L  to  large  values  will  significantly  increase  the  frictional  resistance  in 
the  MHD  duct.  From  Eq.  48  it  can  be  seen  that  increasing  w  decreases  /,  the  total  current 
in  the  duct,  which  increases  r][  in  Eq.  47,  as  the  numerical  results  show. 

The  MHD  propulsive  efficiencies  reported  here  which  approach  those  of  a  propeller 
correspond  to  duct  volumes  greater  than  160  m3,  with  small  surface  area-to- volume  ra¬ 
tios.  These  large-volume  magnetic  inductions  of  6  to  10  tesla  may  be  unattainable  in 
practice. 


MAGNETIC  INDUCTION,  B  =  6  T  »7i 


DUCT  LENGTH,  L  =  10  m 
DUCT  HEIGHT,  h  =  4  m 
DUCT  WIDTH,  w  m  4m 
VEHICLE  SPEED,  V  *  10  m/s 
HYDRODYNAMIC  RESISTANCE  OF 
VEHICLE,  R  -  3,000  V2  N 


0.4  «  CONDUCTIVITY  AT  POLAR  ICE  CAPS 


2.0  4.0 


8.0  10.0  12.0  14.0  16.0  18.0  20.0 


ELECTRICAL  CONDUCTIVITY  OF  SEAWATER,  o  (S/m) 


Fig.  11.  Propulsive  efficiency  versus  electrical  conductivity  of  seawater. 
(Electrical  conductivity  of  normal  seawater  is  approximately  4  S/m.) 
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PROPULSIVE  EFFICIENCY 


Fig.  12.  Propulsive  efficiency  of  square  MHD  duct  versus  duct 
cross-sectional  area. 

Tables  2  to  4  present  the  propulsive  efficiencies  r}X  for  different  length  L  square 
ducts  having  height  h  of  4  m  and  width  w  of  4  m  versus  vehicle  speed  V  in  meters  per 
second.  (Note:  1  m/second  =  1.944  knots.)  The  magnetic  induction  B  across  the  duct  is  2 
tesla  in  Table  2, 6  tesla  in  Table  3,  and  10  tesla  in  Table  4.  The  hydrodynamic  resistance 
of  the  submerged  vehicle  is  3000  V2  newtons  and  thus  requires  a  power  of  3000  V3  watts 
to  propel  the  vehicle  at  a  constant  speed  V.  The  propulsive  efficiency  rj\  increases  signifi¬ 
cantly  with  an  increase  in  magnetic  induction  B.  Bear  in  mind  that  a  magnetic  induction 
of  6  to  10  tesla  across  the  duct  is  a  practical  upper  limit,  imposed  mainly  through  structur¬ 
al  considerations  using  available  superconducting  magnet  technology.  The  efficiency  rj\ 
decreases  significantly  with  decreasing  vehicle  speed  V.  For  the  2-tesla  magnetic  induc¬ 
tion  B  the  efficiency  rj\  tends  io  increase  slightly  as  duct  length  increases  to  40  m.  For 
the  6-tesla  magnetic  induction  the  efficiency  generally  tends  to  increase  slowly  with  duct 
length  L  to  40  m.  With  the  10  tesla  magnetic  induction  B  the  efficiency  rj  \  generally 
tends  to  increase  slowly  with  duct  length  L  to  40  m.  In  these  numerical  calculations  the 
conductivity  of  seawater  a  was  assumed  to  be  4  siemens/m. 
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PROPULSIVE  EFFICIENCY 


Fig.  13.  Propulsive  efficiency  of  square  duct  versus  duct  length. 

Tables  5  and  6  present  the  propulsive  efficiencies  rj  i  for  the  rectangular  duct  of 
width  w  of  4  m  and  height  h  of  4  m  versus  duct  length  L  in  meters  and  vehicle  speed  V  in 
meters  per  second  for  magnetic  inductions  B  of  6  and  10  tesla,  respectively.  The  tables 
show  that  extremely  long  duct  lengths  L  of  over  50  m  generally  degrade  the  propulsive 
efficiency  rjx  with  the  efficiency  decreasing  with  increasing  length.  This  would  be  ex¬ 
pected  since  the  propulsive  efficiencies  decrease  because  of  the  increased  skin  friction  in 
the  duct  with  increased  surface  area. 

Table  7  presents  the  propulsive  efficiency  rj  t  for  the  preceding  MHD  duct  with  a 
length  L  of  10  m  versus  vehicle  speed  V  in  meters  per  second  and  conductivity  of  seawa¬ 
ter  <7  in  siemens/meter.  The  magnetic  induction  across  the  duct  is  6  tesla.  The 
hydrodynamic  resistance  of  the  submerged  vehicle  in  the  seawater  is  3000  V2.  The  pro¬ 
pulsive  efficiency  rj  i  increased  significantly  with  increases  in  fluid  electrical 
conductivity  a  and  decreases  in  vehicle  speed  V.  The  conductivity  of  seawater  can  be 
increased  by  seeding.  However,  it  does  not  appear  that  this  is  practical  for  MHD 
propulsion. 
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Table  8  presents  the  propulsive  efficiency  rj ,  for  the  duct  versus  magnetic  induc¬ 
tion  B  across  the  duct  and  conductivity  o  of  the  fluid  flowing  through  the  duct.  The 
hydrodynamic  resistance  of  the  submerged  vehicle  in  seawater  is  3000  V-  newtons.  The 
propulsive  efficiency  rj ,  increases  significantly  with  increasing  fluid  conductivity  £7 

and  increasing  magnetic  induction  B.  However,  the  conductivity  of  seawater  is  generally 
about  4  siemens  per  meter  and  the  practical  limit  to  B  is  about  10  tesla  with  today’s  pres¬ 
ent  superconducting  magnet  technology. 

Tables  9  through  13  present  the  propulsive  efficiency  rj ,  for  the  square  MHD  pro¬ 
pulsive  duct  versus  vehicle  speed  V  in  meters  per  second  and  magnetic  induction  B  in 
tesla.  The  hydrodynamic  resistance  R  of  the  vehicle  is  1000  V 2  newtons  in  Table  9,  2000 
V2  newtons  in  Table  10,  3000  V2  newtons  in  Table  11,  4000  V2  newtons  in  Table  12, 
and  5000  V 2  in  Table  13.  For  each  table  the  propulsive  efficiency  rj ,  increases  with 

magnetic  induction  B  and  decreases  with  vehicle  speed  V  in  meters  per  second.  These 
table  also  show  that  the  propulsive  efficiency  rj  ,  for  the  MHD  system  generally  increases 

as  the  hydrodynamic  resistance  of  the  submerged  vehicle  decreases. 
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Rg.  15.  Propulsive  efficiency  r|i  versus  duct  height  (duct  width  w  =  2  m). 


Propulsive  efficiency  V|t  versus 
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Fig.  17.  Propulsive  efficiency  Tji  versus  duct  height  (duct  width  w  =  8  m). 


38 


DTRC-90/009 


Fig.  18.  Propulsive  efficiency  rj1  versus 


Table  2.  Propulsive  efficiency  rj ,  x  102  versus  submerged  vehicle  speed  and  length  of  square  duct  at  2  tesla.  (Duct 

height  h  -  4  meters,  duct  width  w  =  4  meters,  magnetic  Induction  in  duct  B  =  2  tesla,  hydrodynamic  resistance  R  =  3000 
V2  newtons,  skin  friction  coefficient  in  duct  Qd  =  0.00262.  Velocity  of  submerged  vehicle  Vin  meters  per  second;  length 
of  duct  L  in  meters.) 
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Table  4.  Propulsive  efficiency  rj ,  x  102  versus  submerged  vehicle  speed  and  length  of  square  duct  at  10  tesla.  (Duct 

height  h  =  4  meters,  duct  width  w  =  4  meters,  magnetic  induction  in  duct  B  =  10  tesla,  hydrodynamic  resistance  R  = 
3000  V2  newtons,  skin  friction  coefficient  in  duct  Cfd  =  0.00262.  Velocity  of  submerged  vehicle  Vin  meters  per  second; 
length  of  duct  L  in  meters.) 
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Table  5.  Propulsive  efficiency?; ,  x  102  versus  submerged  vehicle  speed  and  length  of 

square  duct  at  6  tesla.  (Duct  height  h  =  4  meters,  duct  width  w  =  4  meters,  magnetic  in¬ 
duction  in  duct  B  *  6  tesla,  hydrodynamic  resistance  R  =  3000  V2  newtons,  skin  friction 
coefficient  in  duct  =  0.00262  Velocity  of  submerged  vehicle  Vin  meters  per  second; 
length  of  duct  L  in  meters.) 
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£-25 
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48.2 

40.7 

35.2 

31.0 

27.8 

25.1 

22.9 

21.1 

19.5 

£-50 

55.4 

48.4 

42.9 

38.6 

35.0 

32.1 

29.6 

27.4 

25.6 

24.0 

£-75 

50.0 

44.8 

40.6 

37.1 

34.2 

31.7 

29.5 

27.6 

26.0 

24.5 

£-100 

45.1 

41.0 

37.6 

34.7 

322 

30.1 

28.2 

26.5 

25.1 

23.8 

£-125 

40.9 

37.5 

34.7 

32.2 

30.1 

26.2 

26.6 

25.1 

23.8 

22.6 

£-150 

37.3 

34.5 

32.0 

29.9 

28.0 

26.4 

24.9 

23.6 

22.4 

21.4 

£-175 

34.4 

31.9 

29.7 

27.8 

262 

24.7 

23.4 

222 

21.1 

202 

£-200 

31.8 

29.6 

27.7 

26.0 

24.5 

23.1 

22.0 

20.9 

19.9 

19.0 

£-225 

29.6 

27.6 

25/9 

24.3 

23.0 

21.8 

20.7 

19.7 

18.8 

18.0 

£-250 

27.6 

25.8 

24.3 

22.9 

21.6 

20.5 

19.5 

18.6 

17.8 

17.0 

£-275 

25.9 

24.3 

22.8 

21.6 

20.4 

19.4 

18.4 

17.6 

16.8 

16.1 

£-300 

24.4 

22.9 

21.6 
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16.7 

16.0 

15.3 

£-325 

23.1 

21.7 

20.4 

19.3 

18.3 

17.4 

16.6 

15.9 

152 

14.6 

£-350 

21.9 

20.6 

19.4 

18.4 

17.4 

16.6 

15.8 

15.1 

14.5 

13.9 

£-375 

20.8 

19.6 

18.5 

17.5 

16.6 

15.8 

15.1 

14.5 

13.9 

13.3 

£-400 

19.8 

18.7 

17.6 

16.7 

15.9 

15.1 

14.5 

13.8 

13.3 

12.7 

£-425 

18.9 

17.8 

16.9 

16.0 

152 

14.5 

13.9 

13.3 

12.7 

122 

£-450 

18.1 

17.1 

16.1 

15.3 

14.6 

13.9 

13.3 

12.7 

122 

11.7 

£-475 

17.3 

16.4 

15.5 

14.7 

14.0 

13.4 

12.8 

122 

11.7 

11.3 

£-500 

16.6 

15.7 

14.9 

14.1 

13.5 

12.9 

12.3 

11.8 

11.3 

10.9 

DTRC— 90/009 


Table  6.  Propulsive  efficiency  r\  x  x  10-  versus  submerged  vehicle  speed  and  length  of 

square  duct  at  10  tesla.  (Duct  height  h  =  4  meters,  duct  width  w  =  4  meters,  magnetic  induc¬ 
tion  in  duct  B  =  10  tesla,  hydrodynamic  resistance  R  =  3000  V2  newtons,  skin  friction 
coefficient  in  duct  =  0.00262.  Velocity  of  submerged  vehicle  V  in  meters  per  second, 
length  of  duct  L  in  meters.) 
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40 
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Lm  25 

69.0 

63.0 

58.0 

53.7 

50.0 

46.8 

44Q 

41.5 

39.2 

37.2 

Lm  50 

61.1 

57.8 

54.8 

52.1 

49.6 

47.4 

45.4 

43.5 

41.8 

40.2 

L  *  75 

54.0 

51.7 

49.5 

47.6 

45.8 

44.1 

42.5 

41.1 

39.7 

38.4 

A.  -  100 

48.1 

46.4 

44.7 

432 

41.7 

40.4 

39.1 

38.0 

36.8 

35.8 

L  -  125 

43.4 

41.9 

40.6 

39.3 

38.1 

37.0 

36.0 

35.0 

34.0 

33.1 

L-  150 

39.4 

38.2 

37.1 

36.0 

35.0 

34.1 

33.2 

32.3 

31.5 

30.7 

A-  175 
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34.1 
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32.3 
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30.7 

29.9 
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27.9 

27.2 
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28.0 
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25.5 
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A-  275 

27.1 

26.4 

25.8 

252 

24.6 

24.0 

23.5 

23.0 

22.5 

22.1 

A  ■  300 

25.5 

24.9 

24.3 

23.7 

23.2 

22.7 

22.2 

21.7 

21.3 

20.8 

A  -  325 

24.1 

23.5 

23.0 

22.4 

21.9 

21.5 

21.0 

20.6 

20.2 

19.8 

A  m  350 

22.8 

22.3 

21.8 

21.3 

20.8 

20.4 

19.9 

19.5 

19.1 

18.8 

A-  375 

21.6 

21J2 

20.7 

202 

19.8 

19.4 

19.0 

18.6 

18.2 

17.9 

A-  400 

20.6 

20.1 

19.7 

19.3 

18.9 

18.5 

18.1 

17.7 

17.4 

17.1 

1.-425 

19.7 

19.2 

18.8 

18.4 

18.0 

17.7 

17.3 

17.0 

16.6 

16.3 

A- 450 

18.8 

18.4 

18.0 

17.6 

17.3 

16.9 

16.6 

16.3 

15.9 

15.6 

1.-475 

18.0 

17.6 

17.3 

16.9 

16.6 

162 

15.9 

15.6 

15.3 

15.0 

A  -  500 

17.3 

16.9 

16.6 

162 

15.9 

15.6 

15.3 

15.0 

14.7 

14.4 

DTRC-90/009 


43 


Table  7.  Propulsive  efficiency  T^xlO2  versus  electrical  conductivity  and  submerged  vehicle  speed.  (Duct  height  h  = 

4  meters,  duct  width  w  =  4  meters,  duct  length  L  =  10  meters,  magnetic  induction  in  duct  B  =  6  tesla,  hydrodynamic  resis¬ 
tance  R  =  3000  W8  newtons.  Conductivity  of  seawater  G  in  siemens  per  meter;  speed  of  submerged  vehicle  V  in  meters  per 
second). 
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Table  9.  Propulsive  efficiency  r|,x102  versus  hydrodynamic  resistance  of  submerged  vehicle  (Duct  height  h  =  4  me¬ 
ters,  duct  width  w  =  4  meters,  duct  length  L  =  10  meters,  hydrodynamic  resistance  R  =  1000  V2  newtons.  Velocity  of 
submerged  vehicle  Vin  meters  per  second;  magnetic  induction  Bin  tesla). 
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DISCUSSION  AND  CONCLUSIONS 

A  simple  magnetohydrodynamic  (MHD)  propulsive  configuration  can  be  a  straight 
rectangular  duct  employed  to  produce  the  thrust  for  a  submerged  vehicle.  The  MHD  pro¬ 
pulsive  duct  can  be  within  the  hull  of  the  submerged  vehicle,  in  the  configuration  of  a 
“pod”  attached  by  struts  to  the  main  vehicle,  or  in  an  array  around  the  outside  surface  of 
the  vehicle.  In  the  mathematical  model  discussed  herein,  the  friction  on  the  duct  walls 
due  to  fluid  flow  was  approximately  represented  by  the  ordinary  smooth  pipe  flow  fric¬ 
tional  formulation.  The  vertical  duct  walls  contained  electrodes,  the  horizontal  walls 
contained  superconducting  magnets,  and  the  duct  wall  thickness  was  assumed  to  be  zero. 
Neither  the  power  supplied  to  the  duct’s  superconducting  magnet  configuration  nor  the 
magnet  masses  were  considered  in  the  mathematical  model.  End  effects,  which  degrade 
propulsive  efficiency,  and  bubble  generation,  were  also  omitted. 

The  fundamental  rectangular  duct  MHD  propulsion  system  was  mathematically 
modelled  to  compute  the  propulsive  efficiencies  of  the  system.  The  propulsive  efficiency 
of  the  duct  within  the  hull  of  the  vehicle  was  designated  as  ij  the  propulsive  efficiency 

of  the  duct  as  an  external  pod  was  designated  as  rj  2.  The  propulsive  efficiency 
rj  2  <  rj ,  because  of  the  increased  skin  friction  of  the  system  due  to  the  external  skin  fric¬ 
tion  of  the  pod.  Numerical  solutions  were  presented  for  the  propulsive  efficiencies 
tj  ,  and  jj  2  of  the  MHD  propulsion  duct  in  regard  to  such  variables  as  duct  geometry, 

duct  magnetic  induction,  electrical  conductivity  of  fluid,  submerged  vehicle  speed,  sub¬ 
merged  vehicle  hydrodynamic  drag,  and  skin  friction  fluid  drag  inside  the  duct 

The  propulsive  efficiencies  calculated  from  the  model  in  general  overestimate  the 
efficiencies  of  the  corresponding  real  MHD  propulsion  system.  The  numerical  propul¬ 
sion  efficiency  curves  in  this  report  can  be  used  for  engineering  estimates  of  the 
propulsive  efficiencies  of  real  MHD  propulsion  systems. 

In  this  work,  a  base  MHD  duct  was  chosen  with  width  4  meters,  height  4  meters, 
length  10  meters,  and  magnetic  induction  of.  6  tesla  across  the  duct  propelling  a 
submerged  vehicle  of  3,000  V2  newtons  at  a  vehicle  speed  of  10  meters  per  second.  The 
propulsive  efficiencies  rj ,  and  ij  2  related  to  the  base  MHD  propulsive  duct  propelling 

the  submerged  vehicle  have  the  following  properties: 

•  In  a  6  tesla  field  rj  x  and  rj  2  decrease  from  0.67  and  0.61  at  2.0  meters  per 

second  to  0.24  and  0.22  at  20.0  meters  per  second,  respectively,  a  signifi¬ 
cant  decrease  in  efficiency  with  increasing  vehicle  speed. 

•  At  a  constant  vehicle  speed  of  1 0  meters  per  second  rj ,  and  rj  2  increase 
from  0.07  and  0.06  at  2  tesla  to  0.67  and  0.63  at  14  tesla ,  respectively,  a 
significant  increase  in  efficiency  with  magnetic  induction. 

•  rj ,  and  rj  2  increase  from  0.07  and  0.06  at  an  electrical  conductivity  of  sea¬ 
water  of  0.4  siemens  per  meter  to  0.67  and  0.62  at  conductivity  of  20.0 
siemens  per  meter,  respectively,  a  significant  increase  in  efficiency  with  sea¬ 
water  electrical  conductivity.  If  possible  at  all,  the  conductivity  of  the 
seawater  in  the  duct  could  be  increased  in  practice  only  by  “seeding.”  (The 
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electrical  conductivity  of  seawater  is  on  the  average  approximately  4  sie¬ 
mens  per  meter.) 

•  At  a  constant  vehicle  speed  of  10  meters  per  second  and  a  magnetic  induc¬ 
tion  of  6  tesla,  rj ,  and  rj  2  increase  from  0.14  and  0.13  at  a  duct  cross 

section  of  4  meters2  to  0.62  and  0.53  at  a  cross  section  of  64  meters2 

•  At  6  tesla  rj ,  and  rj  2  increase  from  0.13  and  0.12  at  a  duct  length  of  2.0 

meters  to  0.46  and  0.39  at  duct  length  of  20.0  meters,  respectively,  a  signifi¬ 
cant  increase  in  efficiency  with  moderate  duct  length  at  6  tesla.  Increasing 
the  duct  length  to  200  meters  deceases  rj ,  to  0.30,  as  would  be  expected 

because  of  the  increased  skin  friction  inside  an  extremely  long  duct. 

•  At  vehicle  speed  of  10  meters  per  second,  magnetic  induction  of  6  tesla,  and 
duct  length  10  meters,  rj ,  increases  from  0.24  to  0.43  with  an  increase  in 

the  electrode  height  (i.e.,  duct  height)  from  4  meters  to  32  meters,  a  signifi¬ 
cant  increase  in  efficiency  with  electrode  height. 

•  At  vehicle  speed  of  10  meters  per  second,  magnetic  induction  of  6  tesla, 
duct  height  of  4  meters,  and  duct  length  of  10  meters,  rj ,  increases  from 

0.14  at  a  duct  width  of  1  meter  (i.e.,  distance  between  electrodes)  to  0.60  at 
a  width  of  16  meters,  a  significant  increase  in  efficiency. 

•  At  a  vehicle  speed  of  10  meters  per  second,  rj  x  generally  decreases  with 

increases  in  the  hydrodynamic  resistance  of  the  submerged  vehicle  from 
1000  V2  newtons  to  5000  V2  newtons. 

Thus,  in  summary,  according  to  the  mathematical  model  presented  here  and  the 
range  of  parameters  studied,  the  ideal  MHD  propulsive  duct  efficiencies  increase  in  mag¬ 
nitude  with  increase  in  the  following  factors:  channel  magnetic  induction,  square  duct 
cross  sectional  area  at  moderate  duct  lengths,  fluid  electrical  conductivity,  duct  length, 
distance  between  electrodes,  and  electrode  height. 

The  ideal  MHD  propulsive  efficiencies  decrease  with  vehicle  speed,  and  generally 
decrease  with  increasing  vehicle  hydrodynamic  drag.  Extremely  long  ducts  significantly 
decrease  propulsive  efficiencies  because  of  the  increase  in  internal  duct  friction. 

A  comparison  of  the  ideal  MHD  propulsive  efficiencies  rj  j  and  rj  2  for  a  sub¬ 
merged  vehicle  with  the  corresponding  efficiencies  of  a  conventional  propeller*  (i.e., 
0.85),  generally  shows  the  MHD  efficiencies  to  be  less  under  most  conditions.  The  MHD 
efficiencies  tj ,  and  rj  2  which  approach  that  of  a  propeller  correspond  to  ducts  of  very 

large  volumes  (order  of  magnitude  of  300  meters3  with  small  surface  to  volume  ratios),  or 
require  very  high  duct  field  values  (20  or  25  tesla).  In  large  volumes  of  space  such  as 
MHD  ducts,  magnetic  inductions  of  about  6  to  10  tesla  are  about  the  highest  field 
strengths  possible  with  the  use  of  present  day  superconducting  magnet  construction  tech- 


*  McCarthy.  J.N.,  Head,  Naval  Hydromechanics  Division  at  DTRC  (Code  154),  Private  communication 
(Oct  1989). 
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nology,**  and  producing  these  fields  would  require  a  significant  engineering  effort  in 
magnet  construction. 

MHD  propulsive  ducts  have  no  mechanical  moving  parts  and  consequently  may 
produce  less  acoustic  noise  than  propellers  under  the  same  propulsion  conditions.  Also, 
not  having  a  propeller  removes  the  problem  of  propeller  hull  shaft  seals.  Using  MHD  pro¬ 
pulsion  on  a  submerged  vehicle  to  reduce  the  acoustic  signature  as  compared  with  a 
conventional  propeller-driven  vehicle  is  possibly  a  very  promising  feature,  and  it  should 
be  investigated,  even  though  the  MHD  system  generally  has  a  lower  propulsive  efficiency. 
Therefore,  it  should  be  determined  if  the  MHD  propulsion  technology  is  available  or  can 
be  developed  in  the  future  to  build  a  submerged  vehicle  that  produces  less  acoustic  noise 
then  conventional  submerged  vehicles. 
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**  Superczynski,  MJ.,  Superconductive  Project  Coordinator  at  DTRC,  Private  Communication 
(Apr  1989). 
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APPENDIX  A 

THREE-DIMENSIONAL  EFFECTS 


This  report  presents  a  one-dimensional  control-volume  model  for  an  MHD  duct 
which  is  moving  through  an  unbounded  fluid  and  which  is  accelerating  part  of  that  fluid 
so  as  to  develop  forward  thrust.  In  this  model,  the  velocity  has  only  an  axial  component, 
while  the  velocity  and  pressure  are  uniform  across  each  cross  section  of  the  flow.  Clearly, 
the  flow  upstream  and  downstream  of  the  duct  must  actually  have  three  nonzero  velocity 
components,  not  one;  neither  are  the  pressure  and  velocity  uniform  across  each  cross  sec¬ 
tion.  The  present  model  is  an  approximate  extension  of  the  one-dimensional  model  for  a 
propellor  in  an  open  sea.  The  key  difference  between  the  two  applications  is  that  the 
open-sea  propellor  provides  a  pressure  rise  across  a  disk  of  essentially  zero  thickness  with 
no  other  solid  surfaces,  while  the  MHD  duct  provides  a  pressure  rise  over  the  axial  length 
L  of  the  duct.  The  effect  of  the  finite-length  duct  walls  is  very  similar  to  the  effect  of  a 
cylinder  or  shroud  placed  around  a  propellor.  This  effect  is  treated  in  detail  by 
Kuchemann  and  Weber.25  Here  we  only  summarize  the  key  features  and  their  implica¬ 
tions  for  the  MHD  propulsion  duct. 

The  propellor  alone  in  an  open  sea  creates  a  certain  two-dimensional,  axisymmetric 
flow  with  a  radially  inward  flow  associated  with  the  axial  acceleration.  If  a  cylinder  of 
finite  length  is  placed  around  the  propellor,  it  locally  blocks  the  radial  velocity  compo¬ 
nent.  In  essence,  each  section  of  the  cylinder  acts  as  an  airfoil  in  the  radially  converging 
flow  created  by  the  propellor  alone.  To  be  an  effective  airfoil,  the  flow  must  not  separate 
at  the  leading  edge  of  the  cylinder  and  this  dictates  a  minimum  radius  of  curvature  of  the 
wall  here.  In  other  words,  the  wall  of  the  cylinder  cannot  be  too  thin  at  the  front. 

The  effects  of  the  cylinder  as  an  airfoil  can  be  stated  in  several  equivalent  ways. 
First,  we  can  say  that  the  circulation  associated  with  the  airfoil  section  increases  the  flow 
through  the  propellor.  For  a  given  thrust,  the  increase  in  velocity  is  smaller  with  more 
flow,  so  that  there  is  less  kinetic  energy  lost  in  the  jet  aft  of  the  propellor.  An  ideal  shroud 
always  increases  the  efficiency,  but  an  actual  shroud  adds  skin-friction  drag,  which  to 
some  degree  offsets  the  gain,  depending  on  the  particular  case.  A  second  way  to  describe 
the  effect  of  the  cylinder  is  to  say  that  the  airfoil  produces  a  pressure  on  the  forward,  in¬ 
side  surface  of  the  cylinder  which  is  further  below  the  ambient  pressure  than  the  local 
pressure  would  be  without  the  cylinder.  This  lower  pressure  has  two  effects.  First,  it 
draws  more  fluid  into  the  cylinder,  resulting  in  more  flow  through  the  propellor.  Second, 
the  propellor  still  provides  an  increase  in  pressure,  but  now  the  upstream  pressure  is  fur¬ 
ther  below  ambient  pressure,  so  that  the  exit  pressure  is  not  as  far  above  ambient 
pressure.  As  a  result,  the  jet  aft  of  the  propellor  does  not  accelerate  and  contract  as  much, 
so  again  there  is  less  kinetic  energy  lost. 

The  cylinder  cannot  be  treated  as  an  airfoil  in  the  open-sea  propellor  flow.  The  pres¬ 
ence  of  the  cylinder  changes  the  flow  created  by  the  propellor,  so  that  the  airfoil-cylinder 
and  the  pressure-jump  propellor  disk  must  be  treated  together.  For  axisymmetric 
shrouded  propellors,  Kuchemann  and  Weber25  present  solutions  involving  axisymmetric 
vortex  sheets  which  produce  flows  satisfying  the  boundary  conditions  at  the  shroud  walls 
and  at  the  propellor  disk.  The  thrust  is  now  divided  between  the  propellor  and  shroud  or 
cylinder.  The  propellor  thrust  is  still  equal  to  the  pressure  increase  times  the  disk  area. 
The  pressures  which  are  below  the  ambient  pressure  near  the  leading  edge  of  the  cylinder 
or  shroud  produce  an  additional  forward  thrust  on  the  cylinder.  An  airfoil  alone  would 
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have  no  net  axial  force.  However,  the  propellor  produces  an  abrupt  pressure  jump  on  the 
inside  surface  of  the  airfoil-cylinder.  Thus  the  pressures  on  the  forward  and  aftward  parts 
of  the  airfoil  are  lower  and  higher,  respectively,  than  they  would  be  without  the  propellor, 
and  there  is  a  net  force  on  the  cylinder  or  shroud.  If  the  cylinder  wall  has  zero  thickness 
and  if  the  leading  edge  separation  is  erroneously  ignored,  then  the  pressure  at  the  leading 
edge  is  minus  infinity.  This  infinite  negative  pressure  acting  over  the  zero  area  of  the 
leading  edge  gives  the  finite  forward  thrust  on  the  cylinder. 

The  MHD  equivalent  of  an  open-sea  propellor  would  be  a  uniform  axial  body  force, 
JB,  acting  over  the  rectangular  volume  of  length  L,  but  with  no  duct  walls.  The  stream¬ 
lines  which  would  result  from  such  an  open-sea  body  force  distribution  are  sketched  in 
Fig.  19  The  forward  thrust  on  the  ship  associated  with  this  acceleration  would  be  applied 
to  the  magnet  coils  producing  the  magnetic  field.  When  we  add  the  actual  duct  walls,  they 
locally  block  the  transverse  velocity  components  and  therefore  act  as  airfoils  with  inward 
lift  forces.  There  is  now  an  additional  forward  thrust  on  the  duct  walls,  associated  primar¬ 
ily  with  the  much  lower  pressures  near  the  entrance  associated  with  the  circulation 
produced  by  the  airfoil-duct  walls.  Again  the  duct  walls  would  require  a  finite  thickness 
and  minimum  radius  of  curvature  to  prevent  separation  at  the  entrance.  Indeed,  there 
would  certainly  be  fairings  fore  and  aft  of  the  straight  duct  walls  since  straight  walls 
alone  would  not  represent  the  optimal  airfoil  cross  sections. 

The  two-dimensional,  axisymmetric  analysis  of  Kuchemann  and  Weber25  for  a 
shrouded  propellor  cannot  be  extended  to  the  MHD  duct  in  a  straightforward  fashion. 
First,  the  pressure  jump  across  the  zero-thickness  propellor  is  replaced  by  a  distributed 
body  force.  Second,  the  flow  is  in  general  three-dimensional,  although  some  designs  for 


Fig.  19.  Streamlines  due  to  uniform  body  force  JB  distributed 
over  rectangular  volume  with  no  duct  walls. 
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the  MHD  propulsion  of  submerged  vehicles  might  involve  axisymmetric  flows.  Since  a 
new  two-  or  three-dimensional  analysis  would  by  needed  to  treat  the  effect  of  the  duct 
walls,  it  would  be  inappropriate  to  keep  a  one-dimensional  model  of  the  MHD  body 
force.  One  would  want  to  include  the  effects  of  the  spatially  decaying,  fringing  magnetic 
fields  and  electric  currents  near  the  entrance  and  exit.  We  have  not  presented  such  a  two- 
or  three-dimensional  model  here.  Our  present  purpose  was  to  provide  the  best  predic¬ 
tions  possible  with  a  one-dimensional  model. 

While  we  have  not  determined  the  increase  in  the  flow  through  the  duct  due  to  the 
airfoil  effect  of  the  duct  walls,  we  can  make  a  few  qualitative  statements  based  on  the  re¬ 
sults  of  Kucbmenn  and  Weber.  They  treated  a  propellor  at  the  middle  of  a  long  cylinder 
and  found  that  the  cylinder  increases  the  velocity  at  the  propellor.  The  effect  on  the  ener¬ 
gy  balance  for  a  fixed  ship  velocity  and  total  thrust  is  to  decrease  the  kinetic  energy  lost 
aft  of  the  propellor.  For  most  of  the  cases  considered  here,  the  lost  kinetic  energy  was  rel¬ 
atively  small  compared  to  other  losses,  particularly  Joulean  heating.  Therefore,  even  a 
relatively  large  reduction  on  the  lost  kinetic  energy  may  not  significantly  improve  the 
overall  efficiency  of  the  propulsion  unit.  In  addition,  an  increase  in  the  velocity  inside  the 
duct  increases  the  frictional  losses,  which  vary  as  the  square  of  the  duct  velocity,  so  that 
any  gains  are  offset  to  a  degree  which  depends  on  the  particular  case. 
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